Accurate Quantum Chemical Calculations 
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I. INTRODUCTION 

An important goal of quantum chemical calculations is to provide an under- 
standing of chemical bonding and molecular electronic structure, and over the last 
thirty years this has been largely realized. A second goal, the prediction of energy 
differences to “chemical accuracy” (about 1 kcal/mole), has been much harder to 
attain. First, the computational resources required to achieve such accuracy are 
very large, and second, it is not straightforward to demonstrate that an apparently 
accurate result, in terms of agreement with experiment, does not result from a can- 
cellation of errors. Therefore, in addition to performing very elaborate electronic 
structure calculations, calibration rind the assignment of realistic uncertainties are 
also required. 

Recent advances in electronic structure methodology (1), coupled with the 
power of vector supercomputers, have made it possible to solve a number of elec- 
tronic structure problems exactly using the full configuration interaction (FCI) 
method within a subspace of the complete Hilbert space. These exact results can 
be used to benchmark approximate techniques that are applicable to a wider range 
of chemical and physical problems. They thus provide the necessary calibration 
of existing techniques for electronic structure calculations, and can be used to de- 
termine the origins of any deficiencies. In fact, as we shall show, the calibrations 
indicate that most of the methods for generating many-electron wave functions per- 
form very well, and that a major source of error in even the best calculations arises 
at a more elementary level, in the selection of the atomic expansion basis sets. In 
this review we will discuss the use of FCI wave functions to benchmark simpler 
computational methods, and new approaches to constructing atomic basis sets that 
substantially reduce the basis set errors while keeping the overall computational 
effort manageable. 

In the following Section we will review briefly the methodology of many-electron 
quantum chemistry. In Section III we consider in more detail methods for perform- 
ing FCI calculations, and in Section IV we discuss the application of FCI methods 
to several three-electron problems in molecular physics. In Section V we describe a 
number of benchmark applications of FCI wave functions. In Section VI we discuss 
atomic basis sets and the development of improved methods for handling very large 
basis sets: these are then applied to a number of chemical and spectroscopic prob- 
lems in Section VII, to transition metals in Section VIII, and to problems involving 

2 ORlGsMAL PAGE IS 

OF POOR QUALITY 


PRECEDING PAGE BLANK NOT FILMED 



potential energy surfaces in Section IX. Although the experiences described in these 
sections give considerable grounds for optimism about our general ability to perform 
accurate calculations, there are several problems that have proved less tractable, at 
least with current computer resources, and we discuss these and possible solutions 
in Section X. Our conclusions are given in Section XI. 


II. QUANTUM CHEMICAL METHODOLOGY 

In the clamped-nucleus Born-Oppenheimer approximation, with neglect of rel- 
ativistic effects, the molecular Hamiltonian operator in atomic units takes the form 
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the absence of external fields. The terms in Eq. (1) comprise the electron kinetic 


energy, the nuclear-electron attraction, the electron repulsion, and the nuclear re- 
pulsion, respectively. (Relaxing the assumptions of fixed nuclei and non-relativistic 
motion can be investigated perturbationally, at least for lighter elements, as dis- 
cussed briefly later.) Our goal is to approximate solutions to the time-independent 
Schrodinger wave equation 

HV = EV ( 2 ) 


for this Hamiltonian. Although some limited progress has been made in approaching 
Eq. (2) by analytical methods (2), this is a very difficult procedure that is not yet 
suited to the production of chemical results. Instead, most methods make (implicit 
or explicit) use of basis set expansion techniques: the unknown eigenfunctions of 
Eq. (2) are expressed in terms of a set of n-particle basis functions {$}• Again, 
while it is possible to consider rather exotic functional forms for the $, involving, 
say, interelectronic coordinates, by far the most common approach is to construct 
each $ using a product of molecular orbitals (MOs — one-electron functions) {t/>}: 

<&K = dY[4>i 

1=1 

Here a given function involves an n-fold product of MOs, to which is applied 
some projection operator or operators 0. As electrons are fermions, the solutions to 
Eq. (2) will be antisymmetric to particle interchange, and it is usually convenient to 
incorporate this into the n-particle basis, in which case the $ will be Slater deter- 
minants. The Hamiltonian given in Eq. (1) is also spin-independent and commutes 
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with all operations in the molecular point group, so that projection operators for 
particular spin and spatial symmetries could also appear in 0. The $ obtained in 
this way are generally referred to as configuration state functions (CSFs). 

The molecular orbitals are usually obtained as linear combinations of a one- 
particle basis 

The one-particle basis functions {x} are often referred to as atomic orbitals (AOs). 
The MO coefficients C are obtained by solving an electronic structure problem 
simpler than that of Eq. (2), such as the independent particle (Hartree-Fock) ap- 
proximation, or using a multiconfigurational Hartree-Fock approach ( 1)- This has 
the advantage that these approximations generally provide a rather good estimate 
of the solutions to Eq. (2) — perhaps 99% of the total energy or more thus 
suggesting that an analysis of the many-electron problem and (possibly) computa- 
tional schemes for attacking it can be formulated around them. Lowdin (3) defined 
the “correlation energy” as the difference between the exact energy obtained from 
Eq. (2) and the Hartree-Fock energy, and the term “correlation problem is widely 
used to refer to the problem of computing this energy difference. 

The most obvious use of the n-particle basis {$} in solving Eq. (2) is the linear 
configuration interaction (Cl) expansion 

^ $K C K- (®) 

K 

If the one-particle basis is complete, the use of all possible $ termed complete 
Cl — in Eq. (5) will yield the exact eigenvalue(s) and eigenfunction(s) of Eq. (2). 
Incidentally, as the use of Eq. (5) corresponds to the variational problem of making 
the energy stationary with respect to variations of the coefficients ck, subject to 
normalization of any guess at the ck will yield an upper bound to the true 
energy. In practice, of course, a complete one-particle space is infinite, and so the 
complete Cl problem would also be infinite in dimension. If we choose a finite, 
truncated one-particle space, but approximate ^ as in Eq. (5) using all the possible 
n-particle basis functions, we have a full Cl (FCI) wave function. The FCI wave 
function can be regarded as the exact solution to a Schrodinger equation projected 
onto the finite subspace generated by the truncated one-particle basis. 
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FCI wave functions have many convenient properties, for example, the results 
are independent of any unitary transformation on the MOs used to construct the 
n-particle basis, and as the only approximation made in solving Eq. (2) is the trunca- 
tion of the one-particle basis, we may identify any discrepancies between calculation 
and experiment as arising from this truncation (assuming the Born-Oppenheimer 
approximation and neglect of relativity axe valid). Hence it would appear that 
FCI calculations are an ideal approach for solving the correlation problem. The 
difficulty, of course, is that the factorial dependence of the length of the FCI ex- 
pansion on the number of electrons correlated and the number of MOs used creates 
insuperable computational difficulties for most problems of chemical interest. Even 
by restricting the correlation treatment to a subset of the electrons, by ignoring 
correlation effects involving Is electrons in first-row systems, for instance, the ex- 
pansions will usually be impractically long. To obtain the dissociation energy of 
the N 2 molecule to within 5 kcal/mole of experiment, for example, would require an 
FCI expansion of some 10 14 CSFs, correlating only the ten valence electrons. Thus 
we must develop schemes for truncating both the one-particle and n-particle spaces 
to arrive at computationally feasible wave functions. 

A simple way to implement n-particle space truncation is to use the uncorre- 
lated wave function (which as noted above is a very substantial fraction of the exact 
wave function) to classify terms in the n-particle space. If we consider the Hartree- 
Fock determinant, for example, we can construct all CSFs in the full n-particle 
space by successively exciting one, two... electrons from the occupied Hartree-Fock 
MOs to unoccupied MOs. For cases in which a multiconfigurational zeroth-order 
wave function is required, the same formal classification can be applied. Since only 
singly- and doubly-excited CSFs can interact with the zeroth-order wave function 
via the Hamiltonian in Eq. (1), it is natural to truncate the n-particle expansion 
at this level, at least as a first approximation. We thus obtain single and double 
excitations from Hartree-Fock (denoted SDCI) or its multiconfigurational reference 
analog, multireference Cl (MRCI). 

Truncated Cl methods are only one of the popular approaches to the correlation 
problem. Coupled-cluster (CC) schemes (4) use a different formal ansatz from 
Eq. (5), based on an exponential operator 

= exp(T)^o (6) 

Here T comprises excitation operators that again excite one, two... electrons from 
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MOs occupied in the zeroth-order function tfo to unoccupied MOs. The use of 
the exponential operator guarantees that the energy is additively separable, that 
is, that the energy of two non-interacting systems is the sum of the two separate 
system energies. This holds for Eq. 3) irrespective of what truncation is applied 
to T (say, to only single and double excitations). This size- consistency property is 
not shared by truncated CI. On the other hand, if Eq. (6) is substituted into the 
variation principle the resulting equations for the cluster amplitudes (the weights 
with which the given single, double... excitation terms appear) are too complicated 
for practical solution, and CC methods instead solve a non-linear, non-variational 
equation system for the amplitudes. Hence the CC energy is not a variational upper 

bound. 

Another popular approach to the correlation problem is the use of perturbation 
theory (5). can be taken 35 an unperturbed wave function associated with a 
particular partitioning of the Hamiltonian, perturbed energies and wave functions 
can then be obtained formally by repeatedly applying the perturbation operator 
to 'I'q- Probably the commonest partitioning is the M0ller-Plesset scheme (5), 
which is used where tfo is the closed-shell or (unrestricted) open-shell Hartree-Fock 
determinant. Clearly, the perturbation energies have no upper bound properties, 
but, like the CC results, they are size-consistent. 

A particular advantage of the CI approach, as opposed to coupled-cluster or 
perturbation theoretic techniques, is that it can be readily formulated to handle the 
case in which 'P 0 is multiconfigurational in character. Corresponding multireference 
CC or perturbation theory approaches are much less well developed. As we shall see 
there are many important chemical applications where a multireference approach 
is mandatory, and we shall therefore concentrate mainly on the CI method in this 
article. 

The classification of CSFs into single, double... excitations is straightforward 
and unambiguous for a closed-shell Hartree-Fock reference function. In open-shell 
or multireference cases there are more possibilities for defining these excited CSFs, 
some of which interact with the reference space in the lowest order of perturbation 
theory, and some of which do not. It is very common to exclude the latter exci- 
tations. This is termed restricting the wave function to the first-order interacting 
space (6). 

The choice of reference space for MRCI calculations is a complex problem. 
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First, a multiconfigurational Hartree-Fock (MCSCF) approach must be chosen. 
Common among these are the generalized valence-bond method (GVB) (7) and 
the complete active space SCF (CASSCF) method (8). The latter actually involves 
a full Cl calculation in a subspace of the MO space — the active space. As a con- 
sequence of this full Cl, the number of CSFs can become large, and this can create 
very long Cl expansions if all of the CASSCF CSFs are used as reference CSFs. This 
problem is exacerbated when it becomes necessary to correlate valence electrons in 
the Cl that were excluded from the CASSCF active space. It is very common to 
select reference CSFs, usually by their weight in the CASSCF wave function. Even 
more elaborate than the use of a CASSCF wave function as the reference space 
is the second-order Cl, in which the only restriction on the CSFs is that no more 
than two electrons occupy orbitals empty in the CASSCF wave function. Such ex- 
pansions are usually too long for practical calculations, and they seldom produce 
results different from a CAS reference space MRCI. 

When more than one state of a system is to be investigated, it is possible to per- 
form separate MCSCF calculations, followed by MRCI calculations, on each state. 
This, however, can be a very expensive process, and if transition properties between 
states are desired, such as transition dipole moments for spectroscopic intensities, 
the nonorthogonality between the MCSCF orbitals for the different states creates 
complications. A simple alternative is to perform an MCSCF optimization of a 
single average energy for all states of interest (9). All states are thereby described 
using a common set of MOs. Although these MOs are obviously not optimum for 
any of the states, experience shows this has little effect on the final MRCI results. 

In some cases, even the use of severely truncated reference CSF spaces can give 
rise to an MRCI expansion that is too long for practical calculations. Siegbahn 
has developed an approximate scheme for handling such cases: the weights of the 
singly- and doubly-excited CSFs are estimated perturbationally and frozen, then 
scale factors multiplying the weights of groups of CSFs are optimized in a so-called 
contracted Cl (CCI) calculation (10). Such a calculation requires much less effort 
than the corresponding MRCI calculation, and the results obtained are generally 
good. It may be possible in some cases to perform a large number of CCI calcula- 
tions, say, to characterize a potential energy surface, and to scale the results based 
on MRCI calibration points to improve the accuracy. Another approach to reduc- 
ing the number of variational parameters in the Cl wave function is the internally 
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contracted Cl approach (11). Here single and double excitations are made not from 
individual reference CSFs but from the MCSCF wave function, which is treated as 
a single term. The configurations of the reference CSFs are held at their MCSCF 
values, but this restriction can be relaxed if desired. Internal contraction normally 
produces results similar to the uncontracted MRCI calculation, but at less cost. 

As noted briefly above, truncated Cl expansions do not produce size-consistent 
results. Several approaches of varying sophistication have been developed to address 
this problem. The single-reference Davidson correction (12) 

AE = E CO rr(l ~cl) ( 7 ) 

is a simple perturbation theoretic estimate of the energy effects of higher than 
double excitations. Here E cor r is the correlation energy from single and double 
excitations and c 0 is the coefficient of the Hartree-Fock determinant. Eq. (7) is very 
often extended heuristically to the multireference case (13) using the form 

AE = E corr (l-^c 2 R ). (8) 

R 

Here E corr is the difference between the MRCI energy and the energy obtained 
with the reference CSFs, and c R is the coefficient of reference CSF R in the MRCI 
wave function. We shall use the notation (+Q) to denote the use of the corrections 
given by Eq. (7) or Eq. (8). A problem with such approximations is that as there 
are no useful bounds on the adjusted energy expression some numerical “noise” 
can be introduced across a large set of calculations (say, for a potential energy 
surface), while the MRCI values themselves would be much smoother (14). (This 
“noise” can also be a problem with those MRCI methods that apply selection to 
all CSFs used, rather than just reference CSFs.) A more sophisticated approach 
to the size-consistency problem is to modify the Cl energy expression so that size- 
consistent results are obtained. The earliest such schemes are the various coupled- 
electron pair approximations (CEPA) (15), including the linearized coupled-cluster 
approach of Cizek (16). The single-reference Davidson correction can be viewed as 
a first approximation to the solution of the linearized CC equations. More recently, 
Ahlrichs and co-workers introduced the coupled-pair functional (CPF) (17), which 
has many similarities with the physical reasoning behind the CEPA schemes, and the 
averaged CPF (ACPF) approach (18) which can be extended to the multireference 
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case. 


An advantage of these methods is that it is easier to apply the powerful 
techniques of analytic derivative theory to cases in which a functional is made 
stationary (19). The modified CPF (MCPF) approach of Chong and Langhoff (20) 
is a development of the single-reference CPF method that often gives better results 
than the original when the single- reference description is a poor approximation. 

The wave function that is used to define the reference CSFs is usually the one 
used to optimize the MO coefficients, but the ultimate accuracy of the MOs clearly 
derives from the quality of the AO basis chosen at the outset. In the vast majority 
of current calculations, the AOs are taken as fixed linear combinations of Gaussian- 
type functions centered on the nuclei, but the number of functions used, the range 
of angular momenta, etc, all vary widely among different calculations. We shall 
discuss a number of these aspects in this review, as the one-particle basis set plays 
a fundamental role in determining the accuracy of a given calculation. 

In summary, then, Cl approaches to the correlation problem involve a choice of 
atomic basis set, the optimization of a Hartree-Fock or MCSCF wave function, the 
selection of Cl CSFs (usually by selecting one or more reference CSFs and all single 
and double excitations from them), and the optimization of the Cl wave function. 
There are two potential sources of error in this prescription — the choices of one- 
particle and n-particle spaces. It is these sources of error we shall review in depth 
in the present work. 

III. FULL CONFIGURATION INTERACTION METHODS 

The use of the full Cl (FCI) procedure as an exact test of approximate methods 
has been a key element in improving the accuracy of calculations, and we therefore 
consider the evolution of FCI methodology in some detail. 

In the earliest Cl programs, the Hamiltonian matrix elements were computed 
individually over a list of CSFs that could be chosen more or less arbitrarily (21). 
(This approach is commonly referred to as conventional Cl.) A given Hamiltonian 
matrix element Hkl can be expanded as 

= E E a £ l (pI%) + E E E E b m«(p9 i"), (9) 

V q P 9 r 3 

in terms of MO one- and two-electron integrals (p\h\q) and (pq\rs) and coupling 
coefficients A and B. The latter depend on the occupation and spin coupling of 
the MOs p,q... in the CSFs $ K and $z,. The matrix elements were stored on disk 
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or tape, and after all the elements had been computed the desired eigenvalues and 
eigenvectors were extracted in a subsequent step. Clearly, this approach could be 
used to perform an FCI calculation. However, the factorial increase in length of 
the Cl expansion with number of electrons and orbitals made this possible for only 
very small calculations, such as the special case of three-electron systems. 

The development of the direct Cl approach (22) eliminated the requirement 
for large peripheral storage for the Hamiltonian matrix and greatly reduced the 
time to compute the eigenvalues and eigenvectors, as much redundant work in the 
conventional Cl approach was eliminated. In direct Cl the residual vector 

<TK — T HkL c L ( 10 ) 

L 


is computed from the expression 

-*=£££ WmfrL +£££££ <-(« i«>*. < n > 

L p q L p q r s 

which results from substituting Eq. (9) in Eq. (10). In this way the Hamiltonian 
matrix elements themselves never appear explicitly. However, in the earliest for- 
mulations of direct Cl, each problem (that is, closed-shell single reference, doublet, 
etc) required the development of a new program, as the A and B values in Eq. (11) 
were derived by hand and coded into the program. The first direct Cl program 
for full Cl calculations with more than two electrons was a three-electron FCI pro- 
gram developed by Siegbahn (23). This program was able to handle relatively large 
one-particle basis sets, but three-electron problems are neither common enough nor 
representative enough for such a code to provide many benchmarks. 

An important development in FCI methodology was Handy’s observation (24) 
that if determinants were used as the n-particle basis functions, rather than CSFs, 
the Hamiltonian matrix element formulas could be obtained very simply. Specifi- 
cally, it was possible to create an ordering of the determinants such that for each 
determinant, a list of all other determinants with which it had a non-zero matrix 
element could easily be determined. Further, it was easy to evaluate the coupling 
coefficients A and B in Eq. (11) (the only non-zero values being ±1) and therefore 
to compute the matrix elements. While this greatly expanded the range of FCI cal- 
culations that could be carried out, the algorithm is essentially scalar in nature and 
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therefore rather inefficient on modem supercomputers. Nonetheless, this method 
supplied some important benchmarks, as discussed below. 

Like any general direct Cl method, the shape-driven graphical unitary group 
approach (25) can also be used to perform FCI calculations for some cases, but 
as in Handy’s determinant al scheme matrix elements involving higher than double 
excitations from the reference CSF are processed using scalar algorithms. Like 
the Handy determinantal scheme, this approach has been used to provide several 
benchmarks. 

A major advance in the efficiency of FCI calculations was introduction of a 
factorized direct Cl algorithm by Siegbahn (26). This involves formulating the FCI 
calculation as a series of matrix multiplications: an ideal algorithm for exploiting 
the power of current vector supercomputers. This algorithm is fundamental to our 
present ability to perform FCI benchmarks and we discuss it in detail. We consider 
only the two-electron contribution to a given by Eq. (11), which can be written as 

( 12 ) 

p q r s L 

Siegbahn noted that by inserting a resolution of the identity, B^t can factorized 
explicitly into products of one-electron coupling coefficients: 


p>K L 
n pqrs 


= E4 


JL 


KJ A 

pq ^ rs 


(13) 


This approach had been used earlier to compute the coupling coefficients them- 
selves (27); but if instead the factorization of B is inserted into Eq. (12) the calcu- 
lation of a can be completely vectorized. First, the product of one set of coupling 
coefficients and c is collected in D, 


Di, = £ (14) 

L 

This can be implemented as vector operations; it should be noted that there are very 
few non-zero A values and this sparseness can be exploited in evaluating Eq. (14). 
A matrix multiplication of this intermediate matrix D and a block of the integrals 
is performed, 

E pq = J2(P9\ rs ) D rs- ( 15 ) 


11 



Finally, the intermediate array E is merged with the second set of coupling coeffi- 
cients to form a contribution to <7, 

« = (16) 

J pq 

All steps in the calculation of a are vectorizable, and perform well on vector su- 
percomputers. For large basis sets the matrix multiplication step in Eq. (15) will 
dominate the computational effort. 

This factorized direct Cl scheme still has some difficulties, however. In partic- 
ular, values of the coupling coefficients A are required. While the lists of A values 
are extremely sparse, for very large calculations there will be many non-zero val- 
ues, which must either be precomputed and stored, or computed on the fly. As 
Siegbahn’s intended application was in generating FCI wave functions for use in 
CASSCF optimizations, he stored the coupling coefficients on a disk file. For very 
large calculations disk space becomes the limiting factor, not the CPU power of 
the computer being used. It is not feasible to perform useful FCI benchmark cal- 
culations using this technique in such a form. The modified strategy developed by 
Knowles and Handy (28) overcame this limitation: these authors showed that by 
using determinants instead of CSFs, all non-zero A^ L values (which are ±1) can be 
computed almost trivially on the fly. The Cl vector is much longer when determi- 
nants are used instead of CSFs, but this is not a major handicap when large memory 
computers such as the CRAY-2 are employed. Thus the factorized, determinantal 
FCI scheme allows FCI calculations to be performed for very large expansions and 
hence to be used for benchmarking approximate methods. In essence, the scheme 
arises from a trade-off between memory (and CPU time) and disk storage: such a 
trade-off is not an uncommon feature of programming for modem supercomputers. 

The most time-consuming step in the Knowles and Handy determinantal FCI 
scheme, as noted above, is a series of matrix multiplications, so the code performs 
very efficiently on Cray computers. Calculations as large as 28 000 000 determinants 
have been performed on the CRAY-2 (29); the c and a vectors and some scratch 
arrays must be held in memory, so about 60 million words were required. In fact, 
memory is no longer the resource limitation on current supercomputers, it is disk 
space that becomes critical. The Davidson diagonalization process (30) used to 
obtain the lowest eigenvalue(s) and eigenvector(s) requires the c and a vectors 
from the previous iterations. Therefore, just as storing the coupling coefficients 
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could exhaust the disk storage long before the CPU time became prohibitive in 
Siegbahn’s factorized scheme, so can the storage of the previous c and a vectors in 
the determinantal algorithm. As Davidson noted in the derivation of his iterative 
diagonalization algorithm i 30), it is possible to fold all the previous vectors into one 
vector — in effect, the calculation is restarted with a new trial vector corresponding 
to the best guess available from a fixed number of iterations. While this can slow 
convergence slightly, it limits the disk storage required. In practice, with the CPU 
power of, say, the CRAY-2 and the high degree of vectorization in the FCI scheme, 
as compared to the speed and limits in disk space, it becomes mandatory to use 
folding of the vectors to make the largest calculations feasible. 

Several recent developments should further improve the efficiency of FCI calcu- 
lations. Olsen et al. (31) have derived and implemented a factorized determinantal 
scheme formulated somewhat differently from that of Knowles and Handy. It ap- 
pears that the new scheme should be especially efficient for large basis sets and few 
electrons correlated, which is the situation generally encountered in benchmark cal- 
culations. Another development concerns the form of the n-particle basis. Although 
the use of determinants avoids any storage of coupling coefficients, it generates a 
longer Cl vector than would the use of CSFs and, more importantly, can lead to 
collapse of an excited state to a lower state of a different spin symmetry due to 
numerical rounding in the iterative diagonalization (32). This can be a problem, 
not only in FCI benchmarks, but also in the FCI step in a CASSCF calculation, 
where one may wish to extract many roots in the study of spectroscopic problems. 
Recently, Malmqvist et al. (33) have proposed a method of using CSFs instead 
of determinants, but still avoiding the storage of the coupling coefficients on disk. 
Their approach seems ideal for the CASSCF problem where the number of active 
orbitals is limited, but may not be suitable for the large FCI calculations used for 
benchmarking. 

It is clear that FCI wave function optimization is currently an active area of 
research, and the best method of performing FCI calculations may not yet have 
been achieved. Indeed, in very recent work Knowles (34) has shown how sparseness 
in c and a can be exploited to reduce the computational effort, and Knowles and 
Handy (35) have been able to perform pilot calculations for a 400 000 000 determi- 
nant FCI. Harrison and Zarrabian (36) have modified the factorization used in the 
original Knowles and Handy scheme so as to reduce the dimension of the interme- 
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diate space that appears :n the resolution of the identity: this approach has been 
used to perform FCI calculations with over 50 000 000 determinants. It should be 
feasible to perform calculations with well over 100 000 000 determinants with this 
approach. Thus it may well become possible to perform even larger benchmarks in 
the near future. As we will show below, this is important because for some problems 
there is a coupling of the one- and n-particle spaces, and larger FCI calculations 
will lead to an improved understanding of a larger range of problems. 

IV. FCI APPLICATION CALCULATIONS 

The simplest non-trivial FCI calculation is the three-electron case, for which 
the n-particle problem can be solved exactly for a very large one-particle space. 
Since relativistic effects, which are neglected, are expected to be very small for the 
systems considered in this work, the results should be nearly exact, and hence should 
compare well with experiment. In this section we present three examples where 
three-electron FCI calculations have yielded accurate results of chemical interest. 

The H 2 +H— ►H-{-H 2 exchange reaction has been studied by Siegbahn and 
Liu (37) using the FCI approach in a large one-particle basis set. The potential 
energy surface for this reaction is estimated to be within 1 kcal/mole of the exact 
surface, and has been widely used to evaluate scattering methods. 

FCI calculations have also been used to assess whether H 2 has a positive elec- 
tron affinity (EA) (38). The impetus for this study was a qualitative calculation 
suggesting that H 2 might have a small EA if the “extra” electron attached into 
the 2a g orbital (39). Before asking experimentalists to investigate this possibility, 
it seemed worthwhile to carry out FCI calculations in a large basis set, as this 
requires only a few hours of computer time including convergence tests of the one- 
particle basis. FCI calculations using a large [5s 4 p 3 d If] Gaussian basis set, which 
yields an EA of H that is 0.006 eV smaller than the best value and a D e of H 2 that 
is only 0.03 eV too small, conclusively shows that H 2 does not have a positive EA. 
Thus, it is unlikely that a bound state of H^" can be observed experimentally. 

Another example of a three-electron FCI application is the study of the D 
Rydberg series in the A1 atom (40). The ground state of A1 is P(3s 3 p ). The 
valence occupation 3s 1 3p 2 gives rise to a 2 D state, and it was suggested that one of 
the lowest terms in the Zs 2 nd}( 2 D) Rydberg series might correspond to this occupa- 
tion, although this is contrary to some experimental evidence (41,42). Theoretical 
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studies indicated some mixing of the 3s l 3p 2 configuration into the lower terms of 
the 2 D Rydberg series (43,44). Further, experimental observations by Garton (41) 
suggested that a state lying beyond the ionization limit of the 2 D series is derived 
principally from the 3s 1 3p 2 occupation. Although this problem is amenable to a 
CASSCF/MRCI calculation, the FCI calculations eliminate any concern that there 
is a bias in the treatment of the n-particle space that could lead to an incorrect 
mixing of occupations. FCI calculations in a [7s 6p lid 3/] basis set yield excitation 
energies for the first eight terms in the Rydberg series that are in excellent agree- 
ment with experiment. This agreement as well as direct calculation indicates that 
core- valence correlation has little effect on the term energies, thereby justifying a 
three-electron treatment. Since this FCI calculation solves the n-particle problem 
exactly using a large one-particle basis set, and the term energies agree well with ex- 
periment, the wave functions should reflect the true configurational mixing in these 
states. Although a FCI treatment is independent of the orbital basis, we take the d 
orbitals from the 3s 2 ( x S) state of Al + that corresponds most closely to the Rydberg 
series to simplify the interpretation. Analysis of the FCI wave function indicates 
that the percentage contribution of the 3s 1 3p 2 configuration in the first member of 
the Rydberg series is 24%. The contribution of 3 s 1 3p 2 rapidly decreases so that the 
contributions to the fourth and sixth terms are only 6.1% and <1%, respectively. 
There is 17% of 3s 1 3p 2 occupation not accounted for in the Rydberg series, which 
is concentrated in a state just above the IP limit as suggested by experiment. Thus 
the FCI calculations have definitively resolved the nature of the 3s 1 3p 2 ( 2 D) state. 

V. FCI BENCHMARK CALCULATIONS 

The factorial growth of the FCI treatment with basis set and number of elec- 
trons clearly limits its applicability. However, FCI calculations can be performed 
in moderate-sized (double-zeta plus polarization (DZP) or better) basis sets to cal- 
ibrate approximate methods of including electron correlation. Thus there is con- 
siderable incentive to utilize these FCI methods and the capabilities of modern day 
supercomputers to perform large-scale benchmark calculations. In this section we 
consider several examples of the insight obtained from FCI benchmark calculations. 
In the sections dealing with applications we consider additional FCI benchmark 
calculations in conjunction with specific applications. 

Early FCI benchmark calculations were performed by Handy and co- 
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workers (45,46). Using a double-zeta (DZ) basis set, they considered stretching 
the O-H bond lengths in the H 2 0 molecule to 1.5 and 2.0 times their equilibrium 
values. The FCI results showed that even the restricted Hartree-Fock (RHF) based 
fourth-order many-body perturbation theory (MBPT) approach (47,48), which in- 
cludes the effects of single, double, triple and quadruple excitations, did not accu- 
rately describe the stretching of the bond; the error increased from 0.6 kcal/mole 
at r e to 10.3 kcal/mole with the bonds stretched to twice their equilibrium values. 
Although the MBPT method is rigorously size-consistent and contains the effects 
of higher than double excitations, it does not describe the bond-breaking process 
well, because the RHF reference becomes a poor zeroth-order description of the 
system as the bond is stretched. Size-consistent methods that include double ex- 
citations iteratively — infinite-order methods such as the coupled-cluster approach 
— do better. However, only methods that account for the multireference character 
in the wave function as the bonds are broken, such as the CASSCF/MRCI method, 
provide an accurate description at all bond lengths (49). 

While these earlier studies helped calibrate methods for the approximate treat- 
ment of the correlation problem, there were insufficient studies to draw general 
conclusions about the accuracy of approximate methods. Also, the one-particle 
basis sets were not large enough to eliminate concerns about possible couplings of 
the approximations in the one- and n-particle treatments. Much more extensive 
benchmarking has recently been possible with advances in computers and meth- 
ods (29,32,50-71). Unlike the earlier work, the basis sets are of at least DZP quality 
and a sufficient number of systems have been studied that trends are distinguish- 
able. The FCI studies discussed in the remainder of this work use a version of the 
Knowles and Handy program. As discussed in Section III, this method is based on 
the work of Siegbahn and vectorizes very well. 

The first application of the recent series of benchmarks was the calculation of 
the total valence correlation energy of Ne atom (50). In this work the FCI was 
compared to several single-reference-based approaches. These calculations showed 
that both the +Q correction and the CPF approach gave reasonable estimates for 
the energy lowering of quadruple excitations, and that the accuracy of the different 
treatments varied with basis set. The calculations further showed that higher than 
quadruple excitations were not important in any basis set. It must be remembered 
that Ne is exceptionally well described by a single reference configuration. It is 
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clear from the results of the H 2 O calculations that to obtain the best possible 
agreement with the FCI it will be necessary to go beyond the single-reference-based 
approaches (72). However, the Ne results do suggest that for some problems single- 
reference-based approaches that include an estimate of higher excitations should 
yield reasonably accurate results. 

To evaluate single-reference and multireference Cl approaches to the corre- 
lation problem further, we compare in Table I the FCI 1 A\ - 3 B x separation in 
CH 2 with various truncated Cl results (53). Since the 3 A\ and 3 B\ states are de- 
rived nominally from the 3 P and 5 S states of carbon, respectively, they involve 
different bonding mechanisms that result in a substantial correlation contribution 
to the separation: the SC'F separation is over 14 kcal/mole too large. The error 
of 2.7 kcal/mole at the SDCI level is still much larger than required for chemical 
accuracy (ss 1 kcal/mole). The inclusion of the contribution of unlinked higher ex- 
citations through either the Davidson correction (+Q) or the CPF method reduces 
the error substantially. The origin of the error in the SCF/SDCI treatment is the 
second important configuration, arising from the double excitation 3 aj —* lb\, in 
the 1 A\ state. If the orbitals for the Mi state axe optimized in a two- configuration 
SCF calculation, and correlation is included by performing an MRCI calculation 
based on both these reference configurations, the error is about half that of the 
SDCI+Q or CPF treatments. The error is reduced to only 0.06 kcal/mole if the 
multireference analog of the +Q correction is added. After the 3a \ — ► lb\ excitation, 
the next most important correlation effect is that associated with the C-H bonds. 
If this correlation effect is accounted for in the CASSCF zeroth-order reference for 
a subsequent MRCI calculation, essentially perfect agreement between the MRCI 
and FCI is observed. That is, a well-defined CASSCF /MRCI treatment accounts 
for all of the differential correlation effects. It is interesting to note that adding 
the multireference +Q correction to this MRCI energy results in overestimating the 
effect of higher excitations and the separation becomes smaller than the FCI result. 

We next consider the spectroscopic constants (59) for the ground state of N 2 . 
The values obtained at various levels of correlation treatment with six valence elec- 
trons correlated are summarized in Table II. The SDCI calculation yields a bond 
length that is in good agreement with the FCI, but the error in D e is 0.45 eV, even 
when size-consistency problems are minimized by using the state of N 2 to rep- 
resent two ground state N( 4 S) atoms at infinite separation. Although the addition 
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of quadruple excitations either variationally (SDQCI) or by the +Q, CPF or 
MCPF approximations — further reduces the error in D e , it remains too large for 
chemical accuracy. If both triple and quadruple excitations are included, the spec- 
troscopic constants are all in good agreement with the FCI. However, this level of 
treatment is prohibitively expensive in a large one-particle basis set, and even this 
wave function does not dissociate correctly to ground state atoms, as this requires 
six -fold excitations relative to the SCF configuration at r e . The spectroscopic con- 
stants computed from an MRCI treatment based on a CASSCF wave function are 
in excellent agreement with the FCI. Furthermore, this treatment agrees with the 
FCI for all r values. The addition of the +Q correction does not affect r e or u e , but 
it makes D e too large compared with the FCI. A comparison of the computational 
difficulties encountered in describing the triple bond in N 2 with the single bond 
in CH 2 confirms that electron- dense systems put a much greater demand on the 
computational methodology. 

In addition to the total energy of the system, it is desirable to carry out FCI 
calibration studies of properties such as dipole moments, polarizabilities and elec- 
trostatic forces. For example, in the 0+0H— >0 2 +H reaction (73), the preferred 
approach of the 0 atom is determined by the dipole-quadrupole interaction. At 
long distances this favors a collinear approach to the H atom, whereas for reaction 
to occur the 0 atom must migrate to the 0 end of OH. An accurate description of 
weakly interacting systems such as van der Waals complexes requires a quantitative 
description of dipole-induced dipole or induced-dipole induced-dipole interactions. 
Further, the dipole moment and polarizability functions of a molecule determine its 
infrared and Raman spectral intensities. 

As a first example of an FCI calibration of properties we present in Table III 
a study of the polarizability of F~ (56). As in the previous examples, the SDCI 
treatment contains significant error. The inclusion of an estimate of higher exci- 
tations improves the results; in this case CPF is superior to the +Q correction. 
In the multireference case two different approaches were used. In the first, the 
CASSCF included the 2 p electrons and the 2 p and 3p orbitals in the active space, 
and all CASSCF CSFs were used as references for the Cl, in which the 2s and 2 p 
electrons were correlated. Results obtained in this way are denoted MRCI in Ta- 
ble III. A more elaborate CASSCF calculation, with the 2s and 2 p electrons and 
the 2s, 3s, 2 p and 3p orbitals in the active space, was also performed: the use of all 
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these CASSCF CSFs as references gives results denoted as MRCF in Table III. The 
MRCI and MRCf results do not agree as well with the FCI as do the MRCI+Q 
or MRCI'+Q results. The +Q correction does not overshoot as it did for N 2 and 
CH 2 , in part because of the larger number of electrons correlated here. As noted 
above, when only six electrons are correlated, the CASSCF/MRCI accounts for 
such a large percentage of the correlation energy that the +Q correction overesti- 
mates the remaining correlation. For more than six electrons, or for cases where 
the zeroth-order wave function used is less satisfactory than was the CASSCF for 
N 2 and CH 2 , the +Q correction is a better approximation. This is especially true 
for such quantities as electron affinities that involve large differential correlation 
effects. Thus the +Q correction substantially improves the agreement with the FCI 
for the electron affinity of fluorine (29), even when large CASSCF active spaces and 
MRCI wave functions are employed. 


VI. ANO BASIS SETS 

In the previous section we showed that the CASSCF/MRCI approach yields 
results in excellent agreement with FCI, that is, results that are near the n-particle 
limit. We may therefore expect excellent agreement with experiment when the 
CASSCF/MRCI approach is used in conjunction with extended one-particle basis 
sets. It has become clear (74-76) that, until recently, the basis set requirements 
for achieving the one- par tide limit at the correlated level were commonly underes- 
timated, both in the number of functions required to saturate the space for each 
angular momentum quantum number and in the maximum angular momentum 
required. For the segmented basis sets that are widely used in quantum chem- 
istry (77), improving the basis set normally involves replacing a smaller primitive 
basis set with a larger one. It is then seldom possible to guarantee that the smaller 
basis spans a subspace of the larger set, and it is thus difficult to establish how 
results obtained with different basis sets relate to convergence of the one-particle 
space. Ideally, the possibility of differences in primitive basis sets would be elimi- 
nated by using a single (nearly complete) primitive set, contracted in different ways 
such that the smaller contracted sets are subsets of the larger one. Such an approach 
is best implemented using a general contraction scheme, such as the one proposed 
by Raffenetti (78) for contracting valence orbitals at the SCF level. However, using 
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atomic SCF orbitals to define the contraction is not necessarily suitable for handling 
the correlation problem, and provides no means to contract polarization functions, 
which require large primitive sets for accurate results. Calculations on molecular 
systems have shown (79) that natural orbitals (NO) provide an efficient method 
of truncating the orbital space in correlated treatments. Almlof and Taylor (80) 
have proposed an NO procedure for contracting atomic basis sets suitable for use 
in correlated molecular calculations: this atomic natural orbital (ANO) approach 
is an efficient method for contracting large primitive valence and polarization basis 
sets. It has the advantage that the natural orbital occupation numbers provide a 
criterion for systematically expanding the basis set. 

These ideas axe illustrated for N atom and N 2 in Table IV . As the contraction of 
the (13s 8p 6 d) primitive set is expanded from [4s 3p 2d] to [5s 4p 3d] to [6s 5p 4d], 
the correlation energy systematically converges to that of the uncontracted results. 
The same is true for the (4/) and (2 g) polarization sets. When these ANO sets are 
applied to N 2 , the same systematic convergence of D e is observed. An additional 
advantage of the ANO contractions is that the optimum atomic description they 
provide reduces the problem of basis set superposition error (SE) (81,82), in which 
basis functions on one atom improve the description of another atom, causing a 
spurious lowering of the computed molecular energy. Superposition error is a par- 
ticular problem in calculations that aim at high accuracy, and it is highly desirable 
to minimize its effects. 

In order to treat atomic states with different character equally, the ANOs can 
be averaged to yield a compromise set (83). This is useful for systems with different 
charges, such as F and F~ or transition metals where the radial expectation values 
of the 3d and 4s orbitals are very different for the 3d" 4s 2 , 3d n+1 4s 1 and 3d n+2 
occupations. Generating the average ANOs is analogous to the state averaging used 
to define compromise orbitals suitable for describing molecular states of different 
character. A [5s 4p 3d 2f 1<?] contraction based on the average of F and F has 
an SDCI level EA that agrees with the uncontracted (13s 9 p 6d 4/ 2 g) basis set 
result to within 0.01 eV. This can be compared with a 0.1 eV error for the same size 
basis set that is contracted for F alone, but with the outermost (the most diffuse) 
s and p primitive functions uncontracted. The results are better if the contraction 
is based on F" instead of F, but not as good as using the average ANOs. The loss 
in total energy as a result of averaging is also very small. ANO basis sets averaged 
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for different states should thus supply a more uniform description in cases in which 
there is charge transfer or ionic/covalent mixing. The success for transition metals 
is equally good: the average ANOs yield separations between the lowest states that 
are virtually identical to the uncontracted basis set. It should be noted, however, 
that it may still be necessary to uncontract the most diffuse primitive functions 
and/or add extra diffuse functions to describe properties such as the dipole moments 
and polarizabilities (84,85) that are sensitive to the outer regions of the charge 
density. For example, in Mi atom the polarizability is in excellent agreement with 
the uncontracted basis set. if the outermost s and p functions are uncontacted, but 
half the uncontracted value if these functions axe included in the ANO contraction. 

Dunning (86) has recently suggested that accurate results can be obtained with 
smaller primitive polarization sets than those used in the ANO studies, thereby re- 
ducing the integral evaluation time. This approach is consistent with the ANO 
procedure in that the primitives are optimized at the Cl level for the atoms. Since 
his basis sets contain uncontracted primitive polarization functions and segmented 
valence sets with the outermost functions free, the basis sets are accurate for prop- 
erties as well as the energy. For example, Dunning’s results for OH are in good 
agreement with calculations employing ANO sets. However, these conclusions may 
be system dependent as indicated by the work of Ahlrichs and coworkers on the 
N 2 (87) and F 2 (74) molecules. Despite optimizing the polarization sets at the Cl 
level for the molecules themselves, the energies for these two homonuclear diatomics 
did not converge as quickly with expansion of the primitive sets as with ANOs. We 
conclude that it might be possible to reduce the size of the primitive sets, but this 
requires more study. 

VII. RESULTS FOR SPECTROSCOPIC CONSTANTS 

The FCI benchmark calculations discussed in Section V show that a 
CASSCF/MRCI treatment is capable of accurately reproducing the FCI results, 
at least when six electrons or fewer are correlated. Further, the ANO basis sets 
discussed in Section VI show that it is now possible to contract nearly complete 
primitive sets to manageable size with only a small loss in accuracy. Therefore, a 
six-electron CASSCF/MRCI treatment performed in a large ANO basis set is ex- 
pected to reproduce accurately the FCI result in a complete one-particle basis set, 
and hence should accurately reproduce experiment. For some systems containing 
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eight electrons, the MRCI agrees well with the FCI, while for others that have very 
large correlation effects, such as systems with significant negative ion character, it 
is necessary to include the +Q correction to achieve chemical accuracy. FCI cali- 
bration for systems with more than eight electrons have not been possible, but we 
expect that the +Q correction will improve the results in most cases. In Sections 
VII-IX we illustrate several calculations that have achieved unprecedented accuracy 
by combining FCI benchmarks and ANO basis sets. 

As discussed previously, FCI calculations for CH 2 show that the 
CASSCF/MRCI treatment accurately accounts for the differential correlation con- 
tribution to the CH 2 - 3 Ri separation. In Table V, this level of treatment is 
performed using increasingly accurate ANO basis sets (88). It is interesting to note 
that although the [4s 3p 2d If /3s 2 p Id] basis set contains fewer contracted func- 
tions than the largest segmented basis sets previously applied to this problem (89), 
it produces a superior result for the separation. The largest ANO basis set used 
gives a separation in good agreement with, but smaller than, the T e value deduced 
from a combination of theory and experiment (90). From the convergence of the 
result with expansion of the ANO basis set, it is estimated that the valence limit is 
about 9.05±0.1 kcal/mole. The remaining discrepancy with experiment is probably 
mostly due to core-valence correlation effects. However, as the valence correlation 
treatment is nearly exact, finer effects such as Bom-Oppenheimer breakdown (91) 
and relativity must also be considered. While FCI calculations have shown that a 
very high level of correlat ion treatment is required for an accurate estimate of the 
CV contribution to the separation, theoretical calculations (66) indicate that CV 
correlation will increase the separation by at most 0.35 kcal/mole see later dis- 
cussion. Therefore, it is now possible to achieve an accuracy of considerably better 
than one kcal/mole in the singlet-triplet separation in methylene. 

An analogous study (88) for SiH 2 indicates that the singlet-triplet separation 
can also be accurately computed for this second-row molecule see Table V. In fact, 
the differential valence correlation contribution to T e converges faster with basis set 
expansion than for CH 2 . Note also that the 1 A x state of SiH 2 is better described 
1 v a single reference than is the state of CH 2 . However, it now becomes nec- 
essary to include the dominant relativistic contributions, namely the mass-velocity 
and Darwin terms (92), via first-order perturbation theory (93) or by using an 
effective core potential (94), if chemical accuracy is to be achieved. Once relativis- 
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tic effects and zero-point energy (95) have been accounted for, the theoretical T 0 
value of 20.92 kcal/mole is in excellent agreement with the higher experimental 
value (96). These calculations rule out the alternative value of -18.0 kcal/mole for 
the 1 Ai — 3 B\ separation based on a lower value for the ionization potential of the 
1 A\ state of SiH 2 - Thus in addition to establishing the T e of SiH 2 , these calculations 
also estaolish the IP of the 1 A x state to be 9.15 eV, the higher of two recent exper- 
imental values (see discussion in Ref. 96). We should note that Balasubramanian 
and McLean (97) independently came to the same conclusion. In their calculations 
a more conventional segmented one-particle basis set was used and therefore the 
accuracy of the SiH 2 calculations was not sufficient to definitively determine the 
T e value. However, by treating both CH 2 and SiH 2 to about the same accuracy 
they were able to predict an accurate T t value for SiH 2 using the error in their CH 2 
calculations. The advantage of the ANO basis sets is that they provide a convement 
method of systematically approaching the basis set limit. This reduces the need to 
use an analogous better-known system for calibration. This is advantageous as it 
can be difficult to find such a system for comparison. 

As the calculations on CH 2 and SiH 2 show, theory is now capable of computing 
energy separations between electronic states to better than a few hundred wavenum- 
bers for both first- and second-row systems. Thus, theory has considerable utility 
for determining the ground state when two nearly degenerate states are present, 
which is the case for the Al 2 , Si 2 and N*+ molecules. In the subsequent discussion 
we describe how FCI benchmark calculations combined with ANO basis sets have 
allowed a definitive prediction of the ground state for these molecules. 

Three possible ground states have been suggested for Al 2 based on previous 
theoretical calculations (98-101) and experimental studies (102-105). These are the 
3 £~( 7 r 2 ) and 3 n u (cr* 7 r*) states with two one-electron bonds and the 1 Y,+(cr*) state 
with one two-electron bond. All of these states correlate with the 2 P(3p ) ground 
state of two A1 atoms. The accuracy of the CASSCF/MRCI treatment for AI 2 
was confirmed using FCI calculations in a small basis set (62). In order to confirm 
that the one- and n-particle spaces were not strongly coupled, the FCI calibrations 
were performed in two different one-particle basis sets. For both basis sets the FCI 
and CASSCF/MRCI T e values agreed. CASSCF/MRCI calculations (62), including 
relativistic effects, were then performed in a large ANO basis, showing that the 3 n u 
state of Al 2 lies just 174 cm -1 below the 3 E“ state. The *£+ state was found to 
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lie much higher in energy. While it is straightforward to carry out CASSCF/MRCI 
calculations with six valence electrons correlated, even in very large ANO basis 
sets, it is fax more difficult to include the A1 2s and 2 p electrons at this level of 
correlation treatment. However, the FCI benchmark calculations also showed that 
the CPF state separations were qualitatively correct. Since the CPF method is size- 
consistent, it is expected that increasing the number of electrons correlated will not 
significantly affect its accuracy. It was, therefore, possible to use the CPF approach 
to show that the L-shell and L-M intershell correlation has little effect on the A X 
separation, although it shortens the bond length slightly. Indirect support for the 
theoretical prediction of a 3 n u ground state comes from the failure to observe the 
3 £- <— 3 Ej band system (well known from emission studies) in a jet-cooled beam 
of aluminum clusters (106 ). Further support comes from the work of Cai et al. (107) 
who have observed two absorption bands that they have interpreted in terms of a 
ground state with a vibrational frequency of 284 cm" 1 : this is in good agreement 
with the theoretical prediction of 277 cm 1 for the A 3 II U state (62) and in clear 
disagreement with the known frequency (350.0 cm -1 ) (108) for the A z H g state. 

C 2 , Sh and N^ + are valence isoelectronic, but only for C 2 has there been an 
experimental determination of the ground state (108). Since these molecules con- 
tain multiple bonds, in addition to calibrating the n-particle treatment with the 
FCI approach, we also performed extensive calculations on C 2 to compare with ex- 
periment. Despite the very large correlation energy associated with multiple bonds, 
the CASSCF/MRCI approach was in excellent agreement with the FCI for the 
spectroscopic constants. When this approach was applied to C 2 using an extensive 
ANO basis set, excellent agreement (to within 60 cm x ) was obtained with the 
experimental separations between the A^E^", a 3 n u and b 3 T, g states (71). Similar 
theoretical calculations (71) applied to the Si 2 molecule resulted in a determina- 
tion of the ground state as A 3 E“, and the prediction that the A 3 II U state lies 
440±100 cm -1 higher. The error bars are assigned based on the C 2 calculations, 
the difference between the MRCI and FCI calculations, and the convergence of the 
energy separations with expansion of the ANO basis set. Application of the same 
CASSCF/MRCI treatment to N 2 + (109) yields a X E+ ground state in analogy with 
C 2 . In fact, the a 3 n„ - A X E+ separation is very similar to that found in C 2 , al- 
though the fe 3 Ej state is higher lying in N 2 + - An analysis of both the C 2 and N 2 
wave functions indicates that there is a change in the relative importance of the 
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major CSFs between the CASSCF and MRCI wave functions. This is one reason 
that such high levels of correlation treatment are needed for these systems. The 
FCI calibration calculations were instrumental in showing that the CASSCF/MRCI 
calculations accurately accounted for the effects of valence correlation on the en- 
ergy separations. The ability of theory to accurately treat these systems is a major 
advantage, as spectroscopic determination of the ground states is very difficult. 

In addition to the calculation of accurate separations, the calculation of prop- 
erties is important for an understanding of spectroscopy. We now consider the ac- 
curate calculation of the dipole moment function of the X 2 II ground state of OH to 
determine the strength of the rotational-vibrational bands (Meinel system) (110), 
which are observed in the night sky, in oxygen- supported flames, in the photo- 
sphere of the sun, in OH stars, and in interstellar space. Recently, it has been 
conjectured that the surface-originating glow observed on the Atmosphere Explorer 
Satellites (111) (and possibly the Space Shuttle) is at least partially due to vibra- 
tionally excited OH radic als arising from the collision of the surface with 0( 3 P) 
atoms with a relative translational energy of 5 eV. 

To determine an accurate electric dipole moment function (EDMF) for the X 2 U 
state of OH requires a very high level of correlation treatment, since it is necessary 
to properly account for the 0 — character in the wave function. To calibrate approx- 
imate methods, an FCI dipole moment was computed (61) at five representative r 
values using a [4s 3 p 2d/2s lp] Gaussian basis set. Of the variety of approximate 
methods compared with the FCI, the CASSCF/MRCI treatment reproduced the 
FCI results best, with an error of only 0.001 a 0 in the position of the dipole mo- 
ment maximum. The MRCI spectroscopic constants are also in excellent agreement 
with the FCI. At this level of correlation treatment, it did not make a substantial 
difference whether the dipole moment was evaluated as an expectation value or as 
an energy derivative. Having identified an approximate correlation treatment that 
accurately reproduced the FCI EDMF in a realistic one-particle basis set, this treat- 
ment was then carried out in an extended [6s 5 p 4 d 2/ lp/4s 3p 2<f] ANO basis set. 
This calculation gave a permanent dipole moment that is within 0.01 D of the ex- 
perimental value (112) for u=0, and within 10% of the accurate experimental value 
for the difference in dipole moments between v—0 and v— 1> Rotational-vibrational 
line strengths determine.! from this theoretical EDMF are expected to be more 
accurate than those deduced from a variety of experimental sources. 
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Theory has made a substantial contribution to the determination of accurate 
electronic transition moment functions (TMF) and radiative lifetimes for atomic 
and molecular systems. (See, for example, the recent review articles by Werner and 
Rosmus (113) and Oddershede (114).) However, only recently has it become possi- 
ble to evaluate unambiguously, using FCI benchmarks, the quality of approximate 
correlation methods for determining TMFs. The OH ultraviolet system (A — X) is 
important in many applications, such as combustion diagnostics, and is amenable 
to FCI calculations. Further, there is a large variation in the measured lifetimes, 
which range (with 1 standard deviation) from 625±25 ns for experiments based on 
the Hanle effect (115) to 760±20 ns (116) using the high frequency deflection (HFD) 
technique. The lifetimes measured by laser excitation fluorescence (LEF) range from 
686±14 ns (117) and 693=10 ns (118) to 721±5 ns (119). 

We have studied (63) the convergence of the A — X TM of OH with respect 
to convergence of both the one-particle and n-particle spaces. The FCI calibra- 
tion studies indicated that the state averaged (SA)-CASSCF/MRCI calculations 
reproduced the FCI value to within 0.2%, but that this required including a 6 or- 
bital in the CASSCF active space. Further, a basis set exploration showed that 
the [6s 5 p 4 d 2f lg/4s 3 p 2d] ANO basis set employed in our study is probably 
within 1% of the basis set limit for the transition moment. Thus the computed 
radiative lifetime of 673 ns for the A 2 E + (u'=0, N'= 1) state of OH is expected to 
be a lower bound and accurate to about 2%. This value is in excellent agreement 
with two of the LEF values. The theoretical calculations are sufficiently accurate 
to rule out the somewhat lower value determined by Hanle effect studies, and the 
higher values determined for the t/=0, N '= 1 level by the HFD technique. The 
HFD lifetimes for higher X' values, however, are in relatively good agreement with 
theory — see Fig. 1 where we have plotted the lifetime for t/=0 as a function of 
rotational level. Since the lifetimes are expected to increase monotonically with N ' 
in the regime where the lifetimes are unaffected by predissociation, theory suggests 
that there may be a systematic error in the HFD lifetimes for small N , perhaps 
due to collisional quenching. 

The OH applications discussed here involve transitions between valence states. 
The application of the SA-CASSCF/MRCI approach to valence-Rydberg transi- 
tions uoses additional problems. The requirement of adding diffuse orbitals to the 
one-particle basis set to describe the Rydberg character is well known, so we fo- 
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cus on the requirements for the n-particle treatment ( 68 ). An excellent system 
for such a calibration is the A1H molecule, which can be treated very well as a 
four-electron system, since L-M intershell correlation introduces only a small bond 
contraction (120,121). The FCI study of A1H included the X'Z* and A'Tl valence 
states and the C 1 E + Rydberg state — see Fig. 2. A relatively small basis set was 
used, but one that was capable of describing the Rydberg character of the C state. 
The A 1 n-X 1 S+ transition moment is well described by an SA-CASSCF procedure 
that includes the A1 3s and 3p and H Is orbitals in the CASSCF active space. As 
the C 1 £ + state is derived from the 2 5(3s 2 4s x ) state of Al, the active space must be 
expanded beyond the valence orbitals. When the 4s is added to the active space, all 
MRCI properties except the dipole moment of the C state are in good agreement 
with the FCI. Since the C' 1 E + state is very diffuse, small changes in the shape of 
the Rydberg orbital lead to very large changes in the dipole moment. The CASSCF 
active space must be increased by two <r, two 7 r and one 8 orbital in order to com- 
pute an accurate dipole moment for the C 1 S"*’ state. Thus, the SA-CASSCF/MRCI 
approach is well suited to the study of transitions that involve Rydberg character, 
with only the minor complication that one or two extra orbitals must be added to 
the active space. However, if the Rydberg state dipole moment function is required, 
more substantial augmentation of the active space is needed. 

We now turn to the accurate calculation of dissociation energies (Z) e ), which is 
a more challenging task than either properties or energy separations between states, 
because the calculations must be sufficiently flexible to describe the very different 
molecular and atomic correlation effects equally. Thus only methods that compute 
a very large fraction of the total valence correlation energy can accurately account 
for the differential correlation contribution to D e . We first consider the calculation 
of D e for singly-bonded systems. 

The dissociation energy of NH is of astrophysical interest for the determination 
of its abundance in comets and the sun. It is also of interest in modelling the 
kinetics of rich ammonia flames. However, the experimental determinations since 
1970 have ranged from 3 21 to 3.78 eV and even the recommended values are in 
disagreement (122,123). Quite recently a lower bound of 3.32i0.03 for Do was 
determined using two-photon photolysis of NH 3 (124). The best estimate is obtained 
by combining this lower bound with an upper bound of 3.47 eV, determined from the 
predissociation of the c x n state (125). Thus the D 0 of NH is an order of magnitude 
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more uncertain than that of CH or OH, which leads to unacceptably large errors 
in the astrophysical or kinetic modelling. Since CASSCF/MRCI calculations for 
OH in a large ANO basis set gave a D t value within 0.03 eV of the accurate 
experimental value (108), analogous calculations for the X 3 E state of NH are 
expected to yield results of comparable accuracy. Thus an accurate D e value can 
be computed for NH by using the FCI approach to calibrate the n-particle treatment 
and then applying the same level of theory to CH and OH where accurate D e values 
Eire known experimentally. The FCI calibration calculations (60) showed that the 
CASSCF/MRCI treatment gave equally good D e values for CH, NH, and OH. The 
Do values for CH, NH and OH, obtained using the CASSCF/MRCI approach in 
a large ANO basis set, are summarized in Table VI. A comparison of theory and 
experiment for CH and OH indicates that the calculations underestimate the D 0 
values by about 0.03 eV. Since comparable accuracy is expected for all three systems, 
we can accurately estimate the Do value of NH by adding 0.03 eV to the calculated 
value. The directly computed value cannot be too small, as the ANO basis sets 
have virtually no superposition error, and it is highly unlikely that the error in NH 
is twice that of CH or OH, so we are able to assign an uncertainty of 0.03 eV to our 
recommended value of 3.37 eV. This theoretical prediction is consistent with, but 
just slightly less than, a very recent experimental determination (126). 

The calculation of a D e for a multiply-bonded system or those with many elec- 
trons is even more difficult than for the singly-bonded first-row hydrides such as NH 
or OH. We have carried out a systematic study (127) of the D e values of N 2 , 0 2 , and 
F 2 to evaluate how the errors in the application of the computational methodology 
vary with the degree of multiple bonding. FCI calibration calculations in a realistic 
basis set are only possible for N 2 and 0 2 , and then only with the restriction that 
the correlation treatment is restricted to the 2 p electrons. For both the -V 1 S state 
of N 2 and the X 3 Y,J state of 0 2 , the CASSCF/MRCI treatment correlating the 2 p 
electrons accounts for essentially all of the correlation effects on the spectroscopic 
constants — see Table II. For N 2 this treatment in a large ANO basis set produces 
a D e value that is larger than experiment (108) (see Table VII). Since this basis set 
has virtually no Cl superposition error, this conclusively shows that 2s correlation 
reduces the D e value. Subsequent MRCI calculations including 2s correlation con- 
firm this conclusion, as the MRCI(10) (that is, ten electrons correlated) D e value 
is 0.16 eV less than experiment. The decrease in D e when the 2s electrons are 
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correlated can be explained in terms of an important atomic correlation effect that 
has no analog in the molecular system, namely the 2s — » 3d excitation with a re- 
coupling of the 2p electrons. The fact that the loss of this atomic effect between 
infinite separation and r e could lead to an MRCI(6) value that was too large had 
been suspected (128), but these calculations dramatically demonstrate the effect. 
At the MRCI(IO) level (based on the CASSCF reference space from the six-electron 
calculation — CASSCF(6)) the error in D e is about 4 kcal/mole, or larger than the 
1 kcal/mole desired for chemical accuracy. Recent calculations (129) show that part 
of the error in the MRCI( 10) calculation arises from inadequacies in the n-particle 
treatment: the MRCI reference space must include all possible distributions of the 
2s and 2 p electrons in the 2s and 2p orbitals (note that while the MRCI reference 
space is expanded, the orbitals are taken from the CASSCF(6) calculation), and 
part of the error is due to limitations in the one-particle space. Of the remaining 
one-particle errors, adding more ANOs from the previous primitive set (s through 
g functions) is the most important factor. Further saturation of the s through g 
spaces and higher angular momentum functions ( h and i type functions) are of 
about the same importance as the errors in the n-particle space. Ideally one would 
like to use an FCI calculation to calibrate the errors in the n-particle space, in- 
cluding the N 2s electrons. However, this is not possible currently. Nevertheless, 
the FCI calculations and the systematic improvement of the one-particle basis sets 
obtained using ANO sets give new insight into the requirements for computing the 
D e of N 2 . 

The results for the D e value of 0 2 are in better agreement with experiment 
than for N 2 ; an error of only 1.7 kcal/mole, or only about twice that found for 
the hydrides. For F 2 , FCI calibration calculations are not currently possible, even 
with 2s correlation excluded. However, as in N 2 , the insight obtained from other 
benchmark calculations (29) suggest that to describe the F~F + character in the 
wave function an “extra’" set of 7r orbitals should be added to the CASSCF active 
space. With this addition to the active space, the results are significantly improved: 
even the CASSCF results are now in good agreement with experiment, and the 
MRCI results in a large basis set yield a D e with an error of only 0.5 kcal/mole. 
The agreement is nearly perfect with the addition of the +Q correction. Thus while 
the electron-dense systems are computationally more difficult, relatively accurate 
results can still be obtained using MRCI methods. The accuracy is improved with 
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the inclusion of the +Q estimate for the higher excitations when the number of 
electrons correlated becomes large. Multiple bonds place tighter demands on both 
the one- and n-particle treatments, but the calculations are becoming increasingly 
accurate as the use of larger Cl expansions becomes possible. 

In the remainder of this section we consider problems where the understanding 
of the approximations in the one- and n-particle spaces have allowed important 
chemical problems to be solved. The CN radical is observed in flames, comets and 
stellar atmospheres. The A 2 U - X 2 H + red system is particularly important as a 
nitrogen abundance indicator for red giant stars (130). If the oscillator strength of 
the red system and the ground state dissociation energy of CN are both accurately 
known, the observed emission intensity in a stellar atmosphere can be converted into 
elemental abundances. Unfortunately, there is a large variation in the experimental 
determinations of both quantities. Thus a systematic theoretical study to provide 
reliable values for these quantities was undertaken (131). 

As discussed for N 2 , the computation of accurate bond strengths for multiply 
bonded systems is a difficidt undertaking, requiring both extensive one-particle basis 
sets and a high level of correlation treatment. However, by performing analogous 
calculations for the C 2 , N 2 and NO molecules where the D 0 values are much better 
known, the calculated D„ for CN could be extrapolated by comparing the experi- 
mental and theoretical values for these other multiply-bonded molecules. The s — p 
promotion and hybridization that can occur for C, but not for N or O, could lead 
to a different requirement in the n-particle treatment. The possibility that s — p 
hybridization introduces an error into the computed D e value of CN was excluded 
by computing the D e of CN indirectly from the D t of CN" using the experimen- 
tal EA of C (132) and CN (133). As the two routes agree even though the s - p 
near degeneracy is not present in C“, this effect does not appear to introduce any 
perceptible error in the CN D e value. In this way we were able to arrive at a final 
D 0 value of 7.65±0.06 eV, where the estimated error bars represent 80% confidence 
limits. By increasing the error bars to ±0.10 eV, we believe that the theoretical 
estimate has a 99% probability of encompassing the correct value. 

The calibration of the lifetime of the CN red system was done in two ways: 
the CN violet system, which is well characterized experimentally, was studied, and 
the transitions in the isoelectronic N^ molecule analogous to the red and violet 
system in CN were studied (the Meinel (134) and first negative systems (135)). For 
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the CN violet, the Meinel, and the N| first negative systems the theoretical 
radiative lifetimes agree well with experiment (136-138). Therefore, analogous cal- 
culations are expected to yield an accurate radiative lifetime of the CN red system 
as well. Our computed lifetime of 11.2 /is for the u=0 level of the A 2 II state is 
in excellent agreement with a previous theoretical study (139). Nevertheless, it is 
significantly longer than the value of 8.50±0.45/xs, obtained (140) by extrapolation 
of time decays to zero pressure following photodissociation of C 2 N 2 . It is difficult 
to reconcile this difference, but it is conceivable that there are other decay mecha- 
nisms such as intersystem crossing between the -X’ 2 E + and A 2 II states that shorten 
the experimentally observed lifetime. Other indirect experimental measurements of 
the intensities (141,142), such as measured oscillator strengths, are consistent with 
theory or suggest an even longer lifetime. Note also that the lifetimes deduced from 
analysis of the solar spectrum are consistent with theory if the theoretical Do value 
of 7.65 eV is used in the analysis. Thus the theoretical lifetimes and D 0 value are 
consistent with the solar spectrum model of Sneden and Lambert (130). 

The analysis of comet ary data requires knowing the vibrational transition band 
strengths in the X 2 E+ state of CN. Two very different values for the Einstein 
coefficient (A) of the fundamental 1-0 vibrational band have been reported: one 
based on analysis of cometary data (143) and the other from measurements in 
a King furnace (144). Using the CASSCF/MRCI dipole moment function, the 
computed A x 0 value (145) was in excellent agreement with the value measured in 
the King furnace. The small uncertainty in the computed value suggests that some 
of the assumptions in the model used to analyze the cometary data are in error. 

One of the most important chemiluminescent phenomena is the so-called Lewis- 
Rayleigh afterglow that occurs from the recombination of ground state nitrogen 
atoms (146). Since most of the emission occurs in the first positive bands (1+) 
of N 2 , which cannot be directly populated in the recombination of ground state 
atoms, a precursor state must be involved. Berkowitz, Chupka, and Kistiakowsky 
(BCK) (147) have proposed the A ' 5 E+ state as the precursor, whereas Campbell 
and Thrush (148) invoke instead the A 3 E+ state. The main objection to the BCK 
theory was that the binding energy of the A 1 5 E^" state (estimated (149) at the 
time to be «850 cm -1 ) was not sufficient to maintain an appreciable steady-state 
population. 

Recent theoretical calculations (150) on the A! state potential again il- 
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lustrate that accurate results for the binding energy require both extensive treat- 
ments of electron correlation and a large one-particle basis set. On the basis of 
CASSCF/MRCI calculations using a large ANO basis set, the A! 5 £+ potential 
was found to have an inner well with a depth of about 3450 cm 1 and a substantial 
barrier («500 cm -1 ) to dissociation. Figure 3 shows the B 3 U g , A 3 £+ and 
potential curves and their vibrational levels in the region of interest. The expla- 
nation for most of the experimental observations is apparent from these potentials. 
For example, the barrier in the state yields a quasibound level that allows 

intersystem crossing from the A' state to i/=13 of the B 3 H g state. The observation 
that v' = 13 is populated by the same mechanism as v' < 12, and that emission 
from v' < 13 is similar, but markedly different from v' > 13, even though v' = 13 
is above the energy of the separated atoms, is consistent with the barrier in the 
theoretical potential. At 300 K the maximum intensity originates from t/=ll in 
the B 3 U g state. This arises since the A' state is vibrationally relaxed before m- 
tersystem crossing to the B state, i.e. the lowest vibrational level of the A' state 
most efficiently crosses to u’ = 11 in the B state. At 77 K, an outer van der Waals 
well in the A! state (not shown in the figure), allows for tunnelling to the higher 
vibrational levels of the inner well of the A! state. Since there are fewer collisions 
at this reduced temperature, intersystem crossing to the v' = 12 level of the B 3 U g 
state is more rapid than collisional relaxation, causing the maximum intensity in 
the B 3 B g state emission to increase from v' = 11 to v' = 12. At 4 K the barrier 
in the A' state leads to a cutoff in the emission from v' = 10 — 12. However, the 
A 3 £+ state has no barrier and therefore it populates the B 3 U g state even at 4 K. 
We note that the A 3 £+ and B 3 U g potential curves cross at v = 16 in the A state 
and v'=6 in the B state, thereby giving a maximum in emission for the v'=6 level of 
the B state at 4 K. Therefore, the calculations indicate that the A 3 £+ and A' 5 £+ 
states are both important precursors in populating the B state. The deep well in 
the ■A ,5 £+ state, revealed for the first time by these calculations, eliminates most 
of the objections to the original BCK theory. 

Part of the emission in the Lewis-Rayleigh afterglow involves the Hermann 
infrared system (HIR). This band system consists of a group of unclassified multi- 
headed bands in the region 700-970 nm. Although the system was unassigned, it 
was known (151' that the bands were produced by either a triplet or quintet tran- 
sition and that the system was readily generated from the energy pooling reaction 
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between metastable N 2 (A 3 £+) molecules. As part of our theoretical study (150) of 
the N 2 afterglow, we found that the HIR band positions and intensities corresponded 
exceptionally well to our theoretical values for the C"' 5 II U — A' transition- see 
Table VIII. This, combined with the fact that the energetics are consistent with 
an energy pooling reaction preferentially populating the v'=3 level of the C" 5 II U 
state (152), provided rather convincing evidence for the assignment of the HIR 
system to this transition. This assignment has now been confirmed (153) by a ro- 
tational analysis of the C" 5 Tl u — A' band system produced in emission in a 
supersonic jet discharge. 

VIII. TRANSITION METAL SYSTEMS. 

The accuracy of calculations on transition metal systems has lagged behind 
that of the first and second rows of the periodic table (154). One problem for the 
transition metal compounds is that an SCF wave function is often a much poorer 
representation of the system than for those compounds only composed of first- and 
second-row elements. In fact, for some transition metal systems large CASSCF cal- 
culations are required even for a qualitatively correct description, and a quantitative 
description requires lengthy Cl expansions. In addition to extensive n-particle basis 
requirements, experience has shown that transition metals also require considerably 
larger one-particle basis sets than for the first and second rows. As transition metals 
have such stringent one- and n-particle requirements, FCI benchmark calculations 
and ANO basis sets have given considerable insight into how to improve the accu- 
racy of calculations on transition metal systems. 

A common feature of transition metal calculations is that several atomic asymp- 
totes are involved in the bonding: this mixing is much more prevalent for transition 
metal systems than non-transition metal systems and much harder to describe ac- 
curately. For example, in the transition metal hydrides it has been shown that the 
d n s 2 , d n+1 s 1 and d n+2 asymptotes all contribute to the bonding in the low-lying 
states. The importance of these asymptotes depends on the atomic separations. 
Further, the magnitude of the d - d exchange energy and the relative size of the 
s and d orbitals are both important in determining the bonding in the molecular 
system. 

As a first example of calculations on transition metals, we consider the 
5 D(3d 6 4s 2 ) - 5 F(3d 7 4s 1 ) separation in the Fe atom. FCI benchmark calcula- 
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tions (58) of this separation, in a basis set that includes an / function, show that 
for Fe the single-reference SDCI procedure accounts for most of the differential cor- 
relation energy in the 3d and 4s shells. When this SDCI procedure is carried out 
in a [7s 6 p 4 d 4/ 2 g] basis set and the dominant relativistic effects accounted for, 
the separation is 1.25 eV: this is significantly larger than the experimental sepa- 
ration (155) of 0.875 eV. Adding the 3s and 3 p electrons to the SDCI treatment 
reduces the separation to 1.06 eV. It is highly likely that most of the error in the 
16-electron calculation comes from errors in the n-particle treatment. However, ap- 
proximate methods of accounting for higher excitations cannot be reliably used for 
this problem; the eight electron FCI benchmark calculations showed that the +Q 
correction and the CPF method incorrectly estimate the differential effect of the 
higher excitations. This is due to the different character of 4s and 3d correlation en- 
ergy contributions: for example Cq is reduced much more by 4s — 4s correlation than 
an equivalent amount of 3d — 3d correlation. Thus when two states are compared, 
the +Q correction is reliable only if the state with the larger correlation energy 
also has the smaller c 0 value. The +Q correction may therefore fail to give reliable 
results for transition metal systems when different types of correlation effects are 
included, even though the correction works well in first-row systems when eight or 
so electrons are correlated. This indicates that caution must be used when applying 
approximate treatments that account for higher excitations to single-reference-based 
transition metal calculations. Finally, to compute accurate atomic separations in 
Fe, and presumably for the metal atoms left of it in the first transition row, it will 
be necessary correlate the 3s and 3p electrons. 

As noted above the bonding in transition metal systems frequently involves 
a mixture of several atomic asymptotes. A study (65) of TiH has shown that a 
CASSCF/MRCI treatment agrees well with the FCI for the spectroscopic constants, 
but is much poorer for one-electron properties such as the dipole moment. This 
difference arises because t he bonding in the X 4 4* state of TiH involves a mixture of 
the Ti 3d? 4 s 2 and 3d 3 4s ! atomic occupations. The 3d 3 4s 1 occupation forms a Ti 
4s- H Is bond that is polarized toward the H, thus yielding a large dipole moment 
of Ti + H“ polarity, while in the 3 cP4s 2 occupation, 4s - 4p hybridization occurs, 
with one hybrid orbital polarized toward the H, and the second polarized away from 
the H. The movement of charge toward the hydrogen is balanced by the movement 
away, thus the dipole moment associated with the 3d? 4s 2 asymptote is quite small. 


34 


Since both asymptotes form strong bonds to the hydrogen, but have very different 
dipole moments, the energetics are less sensitive to the mixing than the total dipole 
moment. The FCI studies show that the CASSCF/MRCI method may not yield a 
reliable dipole moment, since the CASSCF generally includes only the nondynanucal 
correlation and the MRCI is not able to fully account for the orbital relaxation 
associated with dynamical correlation. However, the benchmark calculations show 
that the MRCI properties converge to the FCI result if natural orbital iterations 
axe performed, which we denote NO-MRCI. This is different from first-row systems 
where NO iterations were not required; for example, the CASSCF/MRCI properties 
for OH are in excellent agreement with the FCI without NO iterations. Therefore, 
in the study of transition metals at least one NO iteration should be performed to 
test the convergence of properties. 

Since the calculations on TiH showed that the NO-MRCI properties were in 
good agreement with the FCI, it is expected that this level of treatment in a large 
one-particle basis should yield very accurate results. NiH is one of the few transition 
metal systems where the dipole moment (2.4±0.1 D) has been measured (156). At 
the SCF level, the dipole moment of NiH is much too large (about 4 D), because 
the wave function is biased strongly in favor of the 3d 9 4 s asymptote of Ni (157). 
This bias is not fully overcome with low levels of correlation treatment, and at 
the SDCI level the dipole moment (3.64 D) remains much too large. A CASSCF 
treatment allows the Ni -id 8 As 2 and 3d 9 As 1 asymptotes to mix, but because the 
3d 8 4s 2 state lies over one eV too low at this level, the CASSCF mixes too much 
3d 8 4s 2 character into the wave function resulting in a dipole moment that is too 
small (1.74 D). Inclusion of more extensive correlation in the MRCI approach gives 
a better balance between 3d 8 4s 2 and 3d 9 4s 1 and a dipole moment (2.15 D) that is 
closer to experiment. However, the converged NO-MRCI dipole moment of 2.59 D 
is slightly larger than experiment. 

The accurate treatment of weakly bound systems is very difficult, since corre- 
lation is required to describe the dispersion, but superposition errors (SE) can be 
significant at high levels of correlation treatment (82,158). To eliminate the SE, 
very large primitive sets including many sets of polarization functions are required. 
Until recently large valence primitive sets were unavailable for the transition metal 
atoms, making it uncertain whether it was worth studying weakly bound transition 
metal complexes. Partridge (159) has recently optimized large primitive basis sets 
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for the first transition row atoms. When these (20s 12p 9 d) basis sets are supple- 
mented with a diffuse d function to describe the 3d n+1 4s 1 occupation and three p 
functions to describe the 4p orbital, they are of the same quality as the (13s 8 p) 
sets for first-row atoms. These large basis sets have a triple-zeta 4s function, and 
when flexibly contracted and supplemented with extensive polarization sets they 
yield very small superposition errors. Using the ANO procedure these basis sets 
can be contracted to a manageable size. As an example we consider the Ni-H 2 0 
system (85). This study used a (20s 15p 10 d 6/)/[7s 6p 4 d 3/] contraction for Ni, a 
(14s 9p 6d 4/)/[6s 5p 4 d 2 /] contraction for O, and an (8s 6p)/[4s 3p] contraction 
for H. At the SDCI level, the total superposition error is only 0.53 kcal/mole or 
about 10% of the binding energy. The Ni superposition error was only about twice 
that of the water. As discussed in Ref. 85, previous studies of the binding energy 
had obtained spuriously large values as a result of superposition error. Thus with 
the current large primitive basis sets and improved contraction schemes it is now 
possible to study weakly-bound transition metal systems. 

As a further example of accurate calculations for a transition metal system, 
consider the difficulty of determining the ground state of FeH. The laser photode- 
tachment spectroscopic studies of FeH - by Stevens et a l. (160) have been inter- 
preted in terms of a 4 A ground state and a low-lying 6 A state only 0.25 eV higher. 
However, the 4 A — 4 A infrared system, which has been observed in the absorption 
spectra of heated mixtures of iron metal and hydrogen, has not yet been observed 
at low temperatures in matrix studies (161). Although early theoretical calcula- 
tions (162) placed the 6 A state lower, qualitative arguments, based on the r e values 
of the 4 A and 6 A states of FeH and the 5 A ground state of FeH" in relation to the 
width of the peaks in the photodetachment spectra, supported the interpretation of 
a 4 A ground state by Stevens et al . (160). However, more accurate calculations are 
now possible using ANO basis sets and large-scale MRCI calculations (163). For 
example, a CASSCF/MRCI+Q calculation, correlating the eight valence electrons 
in a large [8s Ip 5 d 3/ 2g/As 3 p 2d] ANO Gaussian basis set, predicts the 4 A state 
to be 0.06 eV lower than the 6 A state. (Note that the +Q correction should be 
qualitatively correct since the differential contribution from Cq and the correlation 
energy work in the same direction). The 4 A — 6 A separation is further increased 
to 0.16 eV when a correction for 3s and 3p inner-shell correlation is included based 
on MCPF calculations. Most of the remaining discrepancy with the experimental 


36 



estimate of 0.25 eV for the separation is probably due to underestimating the dif- 
ferential effects of inner-shell correlation. Therefore, these theoretical predictions 
provide strong support for the interpretation of the photodetachment spectra of 
Stevens et al. (160) in which the 6 A state is placed above the 4 A ground state. 
Nevertheless, FeH represents a very difficult case for theory, namely the presence of 
two nearly degenerate states of different multiplicity. Very high levels of correlation 
treatment axe required to quantitatively account for the separation as the correla- 
tion energy is substantially larger for the lower spin state. The failure to observe 
the 4 A - 4 A infrared system in absorption at low temperature is probably due to 
the broadening of this already weak band system by the matrix (160). 

The difficulty in computing an accurate 4 A- 6 A separation in FeH arises mostly 
from the large 3d — 3d exchange energy resulting from the compact nature of the 3d 
orbital. Although the large loss in exchange energy from reversing the electron spin 
is balanced by the greater bonding in the 4 A state of FeH, this latter effect requires 
high levels of correlation treatment to describe. Thus low levels of correlation 
treatment favor the 6 A state, resulting in the incorrect ordering of states. This 
situation can be contrasted with the 1 A\ — 3 Bi separation in CH 2 . More accurate 
results are possible in the latter case, because the differential exchange and bonding 
contributions are much smaller. 

IX. THE CALCULATION OF POTENTIAL ENERGY SURFACES 

In the previous sections we have discussed the application of the most accurate 
theoretical methods for the calculation of spectroscopic constants. Another impor- 
tant area of computational chemistry is dynamics (164). A variety of methods have 
been developed for the determination of scattering cross sections, rate constants and 
product state distributions. These dynamical methods, either classical or quantum 
mechanical, require knowing the form of the potential energy surface (PES). There- 
fore, the accuracy of the predictions of these dynamical methods ultimately depends 
on the accuracy of the PES. It is clear that single-reference-based methods will not 
be able to describe the bond breaking and bond formation that occur during a re- 
action. They may, however, yield accurate surfaces for non-reactive collisions. The 
multireference-based approaches appear to be the only alternative for the accurate 
calculation of general PESs. Given the kind of insight that the FCI calibration has 
given into the n-particle approximations for spectroscopic problems and the high 
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accuracy that the ANO basis sets have shown for these classes of problems, it is 
important to test these advances on the calculation of PESs. Unfortunately it is not 
possible to directly compare a theoretical PES with experimental observations, and 
therefore some of our conclusions will have to wait for confirmation by application 
of dynamical studies to our PESs. 

The first system we consider is the isotope exchange reaction for H 2 and D 2 . 
Experimental studies of this reaction are very difficult to perform since impurities 
allow the formation of H atoms, which then results in a radical exchange reac- 
tion (165). The early theoretical studies were flawed by the assumption that the 
reaction proceeded by a four-center symmetry -for bidden process (166). While the 
later theoretical studies (167) focused on the symmetry-allowed termolecular pro- 
cess, 

2H 2 + D 2 - 2HD + H 2 (17) 

they suffered from limitations in both the one- and n-particle basis sets and con- 
tained an uncertainty of perhaps 10 kcal/mole in the barrier height. However, recent 
calculations (70) have reduced this uncertainty to less than 1 kcal/mole. In that 
study all calculations were performed on six-electron systems in order to minimize 
size-consistency errors in the comparisons, so the H 2 results to be quoted were ob- 
tained from calculations on three separated H 2 molecules. An FCI treatment in 
a DZP basis set shows that the CASSCF/MRCI D e of H 2 is only 0.1 kcal/mole 
smaller than the FCI, while applying the +Q correction to the six-electron calcu- 
lation results in a D e value for each H 2 that is 0.1 kcal/mole too large. There is 
essentially no difference in the saddle-point geometry between the FCI, MRCI, or 
MRCI+Q calculations, and the barrier heights are also in good mutual agreement: 
70.9(MRCI), 70.4(MRCI-i-Q) and 70.4(FCI). The MRCI calculations were then re- 
peated in an ANO basis set that yields a D e value for H 2 at the MRCI+Q level 
(for a system of three separated H 2 molecules) that is only 0.3 kcal/mole smaller 
than experiment. The barrier height at the MRCI and MRCI+Q levels is 66.5 and 
67.4 kcal/mole, respectively. Based on the FCI calibration in the DZP basis set, we 
expect that the true barrier height lies between these values. As the error in the 
barrier height should be smaller than the error for the process of dissociating the 
three H 2 molecules to six H atoms, the barrier should be accurate to better than 
1 kcal/mole, which is comparable to or better than experiment. 

The F+H 2 — ► FH+H reaction has received considerable attention from the- 
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ory (168-174) as it is an excellent system for the calibration of methods: there is 
experimental data for the reaction rate as a function of temperature, information on 
product vibrational energy distribution, and the reaction thresholds are known for 
both H 2 and D 2 (175,176). In Table IX, several different treatments are compared 
to the FCI barrier height and exothermicity (55,168,169); in all of these treatments 
only seven electrons (F 2p and H Is) are correlated. The smallest MRCI treat- 
ment has the F 2 p& and H Is orbitals active in the CASSCF and MRCI (denoted 
MRCI(300), since there are three active orbitals of a x symmetry). This calcula- 
tion yields a barrier height that is 0.68 kcal/mole higher than the FCI value. The 
inclusion of the +Q correction improves the barrier height, but it is now slightly 
too small. The sign of the error in the exothermicity also changes with the addi- 
tion of the +Q correction. Such problems with the MRCI(300) treatment are not 
unexpected since the HF wave function is known to contain significant H + F“ char- 
acter. For accurate results it is therefore necessary to improve the description of 
the electron affinity (EA) of F by expanding the active space to include 2 p* — ► 2 p' n 
correlation. With such a (322) active space, a very large number of CSFs would 
arise in a CASSCF reference MRCI wave function, so it becomes necessary to select 
reference CSFs according to their CASSCF coefficients as described in Section II: 
with this expanded active space the MRCI(322)-(-Q results are in excellent agree- 
ment with the FCI, provided that the threshold for including CASSCF CSFs as 
references is no larger than 0.025. Further expansion of the active space improves 
the results, but the +Q correction is now too large. The ACPF method was also 
used to estimate the effect of higher excitations. The ACPF results are better than 
MRCI(322)(0.025), but not as good as MRCI(322)(0.025)+Q. The differences are 
relatively small, however. It was also found in an ACPF treatment in a large basis 
set that a CSF with a coefficient of less than 0.025 in the CASSCF wave func- 
tion appeared with a much larger coefficient in the ACPF wave function near the 
saddle point. This configuration corresponds to an ionic F - H* contribution, and 
its inclusion (making 13 reference occupations) has little effect in the small basis 
set that does not adequately describe F~. However, it has an important effect in 
the larger basis set, where F — is described much better. In addition, we note that 
singles, doubles, triples and quadruples away from an SCF reference (SDTQCI) has 
a comparable error to the ACPF, but of opposite sign. 

As we noted in the beginning of this section, single-reference-based techniques 
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axe not expected to be quantitative, but they can still be very useful in preliminary 
calculations. In this regard we note that although the CPF method is quite accurate 
for the barrier height and saddle-point geometry, it is significantly poorer than the 
MRCI for the exothermicity. 

In order to compute an accurate barrier height, the basis set is expanded from 
the [4s 3 p ld/2s 1 p] set used for the FCI calibration to a [5s 5 p 3 d 2/ l^/4s 3p 2d] 
ANO set. In this large basis set, the spectroscopic constants for H 2 are in almost 
perfect agreement with experiment. The MRCI(222)+Q treatment of HF, which is 
analogous to the MRCI(322) treatment of F+H 2 , yields an excellent r e , but a D t 
value that is 1.22 kcal/mole (0.05 eV) too small. The Cl superposition error for F is 
0.15 kcal/mole; this is even smaller than that obtained using the large Slater-type 
basis se rom Ref. 170. An analysis of the wave function shows that the saddle 
point resembles F+H 2 more than HF+H. Therefore, the barrier height should be 
accurate to better than 1 kcal/mole. 

The theoretical results for the classical saddle-point and barrier height 
are summarized in Table X. Based on the FCI calibration calculations, the 
MRCI(322)(2p)+Q calculations with seven electrons correlated (i.e. excluding F 2s 
correlation), are expected to reproduce the result of an FCI calculation in a nearly 
complete one-particle basis set. Since F 2s correlation decreases the barrier, this 
MRCI(322)(2p)+Q barrier, when corrected for the Cl superposition error (SE), 
represents an absolute upper bound of 2.52 kcal/mole for the barrier height. The 
ACPF(322)(2p) barrier height is slightly higher than the MRCI(322)(2p)+Q barrier 
height as expected from the FCI calibration calculations. This illustrates that the 
one and n-particle basis sets are not strongly coupled, as both methods of estimating 
the higher excitations agree in two different basis sets. 

The inclusion of F 2s correlation decreases the barrier height, and increases 
the magnitude of the +Q correction. Unlike the seven-electron treatment, the 
ACPF(322) and MRCI(322)+Q results do not agree as well when nine electrons 
are correlated. Unfortunately, it is not possible to calibrate this level of treatment 
using the FCI approach in a realistic one-particle basis set. However, the ACPF 
wave function shows that there is an F H^" contribution to the wave function that is 
not accounted for by simply adding the extra ir orbitals to the active space. As noted 
previously, it is necessary to add an additional CSF to the reference list to bring the 
MRCI and ACPF barrier heights into agreement. The additional CSF significantly 
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changes the ACPF barrier, but does not affect the MRCI barrier and increases the 
MRCI+Q barrier only slightly. The agreement of the two methods gives strong 
support for a conservative estimate of 2.00 kcal/mole for the barrier height, after 
correction for SE. A more realistic value for the barrier height should account for 
basis set incompleteness and the underestimation of the effects of higher excitations 
by the +Q correction: we assume that the +Q estimate is too small, since the 
difference between the size-consistent ACPF(322)(2p) and ACPF(322) treatments 
is larger than the difference between the MRCI(322)(2p)+Q and MRCI(322)+Q 
treatments. Our best empirical estimate of 1.60 kcal/mole is obtained by omitting 
the SE correction, adding instead 0.1 kcal/mole for basis set incompleteness and a 
further correction of 0.1 kcal/mole for an understimation of the + Q reduction to the 
classical barrier height. Thus based solely on estimates from ab initio calculation, 
the barrier height should be between 1.60 and 2.00 kcal/mole. 

Our estimate is very similar to that made by Truhlar and co-workers (171) us- 
ing the scaled external correlation (SEC) method, but different from that suggested 
by Schaefer (174) (i.e. a barrier height greater than 2.35 kcal/mole) or the Monte 
Carlo calculations (172) that yield a barrier of 4.5±0.6 kcal/mole. The lower bound 
for the barrier suggested by Schaefer (174) appears to be supported by SDTQCI 
results (173) (correlating only 7 electrons), that yield a barrier of 2.88 kcal/mole. 
Further, no differential effect of F 2s correlation is observed at the SDTQCI level in a 
small basis set. However, the basis sets used in the SDTQCI calculations are far from 
the one-particle limit. If an estimate for errors in the one-particle basis is included, 
the SDTQCI results are found to be in good agreement with the MRCI(322)(2p)+Q 
and ACPF(322)(2p) methods, as expected based on the FCI calibrations. There- 
fore, the only difference between the MRCI (or ACPF) and SDTQCI is the effect 
of F 2s correlation on the barrier. In Table XI we give the results of FCI bench- 
mark calculations correlating nine electrons in a DZ+pp (double-zeta plus diffuse 
p on fluorine) basis. As we discussed above, for realistic benchmarking a larger 
basis set is required. This is clear from the fact that at the seven-electron level in 
the larger basis set the barrier height is MRCI>ACPF>MRCI+Q>FCI>SDTQCI, 
while in the smaller basis set it is SDTQCI>MRCI>FCI>ACPF>MRCI-(-Q- These 
calculations are only qualitative, but they do suggest that the SDTQCI calcula- 
tion underestimates the differential effect of the F 2s correlation. Thus, even a 
single-reference-based technique that correctly includes both triple and quadruple 
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excitations is inferior to the CASSCF /MRCI approach. We also feel that the defi- 
ciencies of the Monte Carlo method are related to the use of a single reference to fix 
the nodes: the ab initio results clearly show that the wave function near the saddle 
point has a much greater multireference character than the wave function for either 
the products or reactants. 

While the MRCI(322)+Q (or ACPF(322)) calculations in the ANO basis set 
are more reliable than any previous results, considerable computer time would be 
required to compute a global surface at this level. The barrier height was therefore 
investigated using the contracted Cl (CCI) approach. In the same large ANO basis 
set, the CCI+Q barrier is 0.4 kcal/mole higher than the corresponding MRCI+Q 
value. Modification of the basis set was also investigated at the CCI level: / po- 
larization functions on H were found to lower the barrier by only 0.07 kcal/mole, 
while eliminating the g function on F increased the barrier by 0.12 kcal/mole. These 
observations are consistent with the contention that the basis set is effectively com- 
plete. The CCI calculation in this basis set is sufficiently inexpensive that much 
larger regions of the PES can be investigated. Of course, given the difference be- 
tween the MRCI and CCI barrier heights, some account would have to be taken 
of the errors in CCI treatment; this might involve adjusting the parameters in the 
fitted potential based on the MRCI(322)+Q calculation or on information deduced 
from experiment. 

While the CCI+Q PES scaled using the MRCI(322)+Q results should be ac- 
curate, direct comparison with experiment is difficult. To facilitate comparison we 
have employed canonical variational transition state theory (177) at the classical 
and adiabatic barrier using the CCI+Q potential for both F+H 2 and F+D 2 . These 
calculations account for the zero-point energy and include a tunnelling correction. 
The results of these calculations are summarized in Table XII. As expected, the 
zero-point and tunnelling corrections are different for H 2 and D 2 . At the classical 
saddle point, the barrier heights for H 2 and D 2 differ by 0.2 kcal/mole, whereas 
at the adiabatic saddle point the barriers are essentially the same. The observa- 
tion (175,176) of nearly identical thresholds for H 2 and D 2 also provides strong 
support for using the adiabatic barrier. In order to bring the computed threshold 
into agreement with experiment (175,176), we must lower the CCI+Q classical bar- 
rier by 0.7— 0.8 kcal/mole. This produces a barrier height of 1.3— 1.4 kcal/mole, 
or after accounting for the errors associated with various approximations, a barrier 
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height of 1.0- 1.5 kcal/mole. This is in good agreement with the estimate made 
directly from the various multireference treatments, and also with that deduced in 
recent calculations by Truhlar and co-workers (171), although it disagrees with the 
value inferred by Schaefer (174) from most previous calculations. 

The previous two applications consider dynamical problems where only the 
ground state potential energy curve is of interest. However, this is not always the 
case — many situations involve a curve crossing. While most scattering formalisms 
are developed in a diabatic representation, a theoretical PES is computed in the adi- 
abatic representation. Hence when curve crossings (or more complicated phenom- 
ena for polyatomic systems) occur both potentials must be accurately represented 
in the crossing region, and nonadiabatic coupling matrix elements (NACMEs) will 
be required to define the unitary transformation between the diabatic and adia- 
batic representations. Until recently, NACMEs were computed either using finite 
difference methods (178) or via approximations to avoid computing matrix ele- 
ments between nonorthogonal wave functions (179,180). However, Lengsfield, Saxe, 
and Yarkony (181,182) have recently developed an efficient method of evaluating 
NACMEs based on state- averaged MCSCF wave functions and analytic derivative 
methods. This should provide NACMEs of similar overall accuracy to that obtained 
for the adiabatic potentials. 

To accurately describe curve crossings requires, in addition to the NACMEs, 
an equivalent treatment of both states involved in the crossing. In curve crossings 
where the molecular orbitals for the two states are similar, such as interactions be- 
tween valence states derived from diff mt asymptotic limits, the CASSCF/MRCI 
approach would be expected to describe both potentials accurately irrespective of 
which state is used for the orbital optimization. However, when the character of 
the two states is very different, such as valence/Rydberg mixing (68) or interac- 
tion between states derived from ionic and covalent limits (69), it is more difficult 
to achieve equivalent accuracy for the lowest adiabatic state on either side of the 
crossing point. This is commonly the case for charge-exchange reactions, such as 
N++N 2 -+N+N+, or chemi-ionization processes such as M + X — * M + + X , where 
the optimal molecular orbitals for the ionic and neutral solutions differ greatly. To 
gain additional insight into the computational requirements for describing the po- 
tentials in the region of curve crossings, we have studied (69) the Li + F — ► Li + +F 
chemi-ionization process using the FCI approach. In LiF, the lowest adiabatic state 
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at short r values, namely "he ionic A 1 S ^ state, dissociates adiabatically to neutral 
ground-state atoms. There is an avoided crossing at the point where the energy 
difference between the F electron affinity (EA) and the Li ionization potential ex- 
actly balances the 1/r electrostatic stabilization. Since the CASSCF description 
of F“ is poor (29), the CASSCF estimate for the bond distance at the crossing 
point is unrealistically small. When orbitals from the ground-state CASSCF wave 
function are used to construct an MRCI wave function, the CASSCF description 
of the crossing point will compromise the accuracy of the MRCI description. This 
problem is not easily resolved by expanding the CASSCF active space, as very 
large active spaces are required to obtain a good description of the atomic electron 
affinities. However, by performing instead an SA-CASSCF calculation, in which 
the orbitals are optimized for the average of the two lowest 1 E"*’ states in LiF (the 
ionic and neutral states), the orbital bias is eliminated and the MRCI treatment 
is in excellent agreement with the FCI. It is also important to note that this av- 
eraging does not significantly degrade the description of the system near r e . Thus 
state averaging appears to be an excellent method of achieving equal accuracy for 
two potential curves in a curve crossing region, and should also perform well for 
polyatomic systems. The utility of state averaging as a means of obtaining a good 
compromise set of orthogonal molecular orbitals for use in an MRCI wave func- 
tion has also been found to be an excellent route to computing accurate electronic 
transition moments (63), as noted above. Thus there is much in common between 
methods that account accurately for differential correlation effects on a PES and 
those that yield accurate spectroscopic constants and molecular properties. 

Application of the SA-CASSCF/MRCI approach to the LiF curve crossing in a 
(13s 9p 6 d 4/) /[4s 3p 2d If] F and (14s 7p 4d)/[5s 4 p 2d] Li basis set yields smooth 
dipole moment and potential energy curves, as expected. These calculations were 
able to demonstrate that t he previous failure of an analysis of the curve crossing was 
due to poor adiabatic potentials, not the use of the qualitative Rittner model (180). 
The NACMEs deduced from the SA-CASSCF/MRCI calculation are in good agree- 
ment with other determinations (183). The limitations in these calculations arise 
from the difficulty in computing the EA of F. As we discuss in the next section, 
this is still a very challenging task. As our ability to compute electron affinities im- 
proves, the accuracy of such ionic-covalent curve crossings will also improve, as the 
FCI benchmarks have shown that the SA-CASSCF/MRCI calculations are correctly 
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describing this region of the potential. 

The applications considered in this section illustrate that current methods are 
yielding more accurate potential energy surfaces, although it is often too expensive 
to compute entire PESs using very extensive MRCI calculations. However, with 
advances in surface fitting techniques it may be possible to merge accurate calibra- 
tion calculations at critical points on the PES with global calculations using less 
computationally expensive methods. For example, we have used the CCI method in 
a smaller basis set to map out the F + H 2 potential (168). This more approximate 
level of treatment agrees with the best calculations to within about 1 kcal/mole and 
is therefore expected to be semi-quantitatively correct. Scaling of these results can 
thus be expected to yield accurate global surfaces. 

X. SYSTEMS WITH A STRONG COUPLING OF THE ONE- AND n-PARTICLE 
SPACES 

In the previous sections we have illustrated how the accuracy of ab initio cal- 
culations has been greatly enhanced by recent developments. However, there are 
a few systems where even these new advances have not fully answered the ques- 
tion as to which is the appropriate quantum chemical treatment. These problems 
involve correlation of core electrons and computation of electron affinities and we 
first consider core correlation. 

As noted above, a high level of correlation treatment is required to fully account 
for the effect of 3s and 3 p correlation on the 5 D - 5 F separation in Fe atom (58). 
It is impossible to study inner-shell correlation effects for transition metal systems 
using FCI benchmark calculations, and furthermore errors associated with other 
approximations are generally as large. However, the calculations for the 1 A\ — B\ 
separation in CH 2 (88) indicate that the treatment of the valence correlation is 
nearly exact and one of the largest remaining errors is the neglect of correlating the 
C Is electrons. In addition, as correlating the core in C adds only two extra electrons 
it is amenable to FCI benchmark calculations. Thus we have investigated (66) the 
effect of core correlation on the 1 A\ — * B\ separation in CH 2 , as well as the C 3 P — 5 S 
separation. Correlating the core not only places more stringent requirements on the 
n-particle treatment, it also increases the demands on the one-particle basis set. 

The study of the error associated with the ANO contraction procedure on the 
3 P — 5 5 separation in the C atom shows that while the ANO procedure works 
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very well for valence correlation, the convergence of the basis set with contraction 
level for the core-core (CC) and core-valence (CV) correlation is very slow. This is 
related to the fact that the change in c 0 for a given correlation energy lowering, and 
hence the size of the natural orbital occupation numbers associated with the core 
correlation, is smaller than those associated with valence correlation energy. Thus, 
when both core and valence correlation are included, it becomes impossible to ordc 
the NOs based on their contribution to the correlation energy using the occupation 
numbers, and therefore it is difficult to use the ANO procedure effectively. This 
is analogous to the origin of the failure of the +Q correction in describing higher 
excitations in Fe atom. The contraction problem is reduced somewhat when the 
nearly constant CC correlation is deleted and only the CV correlation is included, 
but the convergence with contraction level is still much slower than when only 
valence correlation is included. Note that elimination of the CC correlation also 
reduces the requirements on the primitive basis set. 

In addition to problems associated with the one-particle basis set, the FCI 
benchmarks show that very high levels of correlation treatment are required to 
correctly balance the CV (or CC) and valence correlation. This usually involves 
increasing the active space such that a larger fraction of the valence correlation 
energy is included in the zeroth-order wave function. In fact, for an MRCI calcu- 
lation to account for all of the CV effect, the zeroth-order wave function has to 
account for essentially all of the valence correlation. This leads to a large expansion 
in the valence part of the problem, which when coupled with the larger basis set 
requirements and additional electrons being correlated, leads to a much larger com- 
putational problem than the study of valence correlation only. Failing to expand the 
valence treatment adequately can lead to an overestimation of the differential CV 
contribution even though the total CV correlation energy may be underestimated. 

The effect of CV correlation on the 1 A\ — 3 B\ separation in CH 2 was estimated 
by using a totally uncontracted C s and p basis set to avoid the contraction problem. 
Using the C 3 P — 5 S separation as a calibration, it is likely that the largest possible 
valence active space that could be used in the CH 2 calculations would lead to a 
slight overestimation of the CV effect. However, even this calculation resulted in 
a Cl expansion of about 1.5 million CSFs. The benchmark studies (66) of CV 
correlation indicate that it increases the 1 Ai — 3 i?i separation in CH 2 by about 
0.35 kcal/mole — see Table V. Unfortunately, such CV calculations are too large 
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to be routinely carried out for more complex systems. Thus there is considerable 
motivation for developing less expensive ways of including this effect, such as the 
CV operator approach being developed by Meyer and co-workers (184). 

The nuclear hyperfine interaction, which depends on the spin density at the 
nucleus, is very sensitive to CC and CV correlation. To determine even qualitatively 
correct results for the nitrogen atom hyperfine coupling constant requires correlat- 
ing the Is orbital (185). Very high levels of correlation treatment as well as large 
flexibly contracted basis sets are required for quantitative results. Further there is 
a strong coupling of the one- and n-particle basis sets, which makes it difficult to 
perform definitive FCI calibration calculations (67). To achieve quantitative agree- 
ment between the MRCI and FCI hyperfine constants for N atom, it was necessary 
to include the 2s, 2 p, 3s, 3p and 3d orbitals in the CASSCF active space. This 
is consistent with the calculations to determine the C 3 P - 5 5 separation, which 
showed that the zeroth-order reference must account for almost all of the valence 
correlation energy in order for the MRCI to accurately account for CV and CC 
effects. Also, the basis set must include several / functions to be within 5% of 
experiment. In spite of the strong coupling between the basis set and correlation 
treatment, the FCI calculations demonstrated that diffuse s functions contribute 
significantly to the spin density, because of an important class of configurations in 
volving 2s — + 3s promotion with a spin recoupling. Thus, the study of the isotropic 
hyperfine constant of N atom indicates that quantitative agreement with experiment 
requires very large one-particle basis sets (flexibly contracted in the core and aug- 
mented with diffuse and polarization functions), as well as an extensive treatment 
of electron correlation (i.e. one that correlates the core and accounts for essentially 
all of the valence correlation energy). 

The last application that we discuss that has a strong coupling between the 
one- and n-particle spaces is the calculation of the electron affinities (EAs) of atomic 
oxygen (52) and fluorine (29). We have noted earlier for OH and F+H 2 that extra 
orbitals have to be added to the CASSCF active space to describe the negative 
ion character. The determination of accurate EAs is difficult because the negative 
ions have more correlation energy than the neutrals and the correlation effects are 
different. The one-particle requirements are also greater for the negative ions in 
that the basis sets must be augmented with diffuse functions. For example, in 
0/0“ if only the 2 p electrons are correlated, improving the level of treatment from 
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an SCF/SDCI to a CASSCF/MRCI that includes the 2 p and 2 p' orbitals in the 
active space increases the EA from 1.06 to 1.26 eV. Further expansion of the active 
space, to include the 3 d orbital, increases the EA by only 10% of the effect of 
adding the 2p'. When the 2s and 2 p electrons are correlated, the change in the 
EA is even larger when tlie 2 p' orbital is added to the active space, from 0.99 at 
the SCF/SDCI level to 1.28 eV at the CASSCF/MRCI level. In fact, at the SDCI 
level the inclusion of 2s correlation decreases the EA, while at the CASSCF/MRCI 
level 2s correlation increases it. It is highly likely that the effect of adding the 
3d orbital to the active space will increase when both 2s and 2 p correlation are 
included. Thus, there is a strong coupling between the importance of 2s correlation 
and the level of correlation treatment. Less extensive studies on F/F~ suggest that 
the +Q correction helps in accounting for this 2s effect, but large reference spaces 
are still required. The dimension of the n-particle problem is further increased by 
the fact that tight thresholds must be employed in selecting the reference space 
based on the CASSCF wave functions to avoid significant discrepancies with the 
unselected CAS reference space MRCI. Thus the study of atomic EAs indicate that 
as molecular systems become increasingly ionic, larger theoretical treatments are 
required to obtain equivalent accuracy for the neutral and ion. 

For the three applications discussed in this section, the FCI benchmark calcu- 
lation as well as the CASSCF /MRCI calculations in the large basis sets have given 
insight into the reason why such extensive calculations are required. Unfortunately, 
this level of treatment is not routinely possible. It is possible that a CV operator 
approach (184) will allow the inclusion of this effect without greatly expanding the 
calculation. Since the +Q correction appears to improve the results for the systems 
considered in this section, it is also possible that improved methods of estimating 
the higher excitations will eliminate the need to expand the zeroth-order valence 
treatment over that which is required for only the valence correlation. 


XI. CONCLUSIONS 

The cumulative experience derived from FCI benchmark calculations indicates 
that the more sophisticated n-particle space treatments in use, particularly the 
CASSCF/MRCI method, generally provide a close approximation to the solution of 
the correlation problem. In cases in which bonds are not broken, or excited states 
are not considered, in light atom molecules, single-reference-based treatments of 
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exact or approximate coupled-cluster type also generate high-quality wave func- 
tions. With current supercomputers it is perfectly feasible to consider CC or MRCI 
expansions of more than 10 6 terms, which opens up large areas of chemistry and 
spectroscopy to accurate investigation. Indeed, the main factor that determines the 
accuracy of such treatments is the completeness of the one-particle space. Although 
there has been justifiable pessimism in the past about the slow convergence towards 
basis set completeness with respect to both radial saturation and angular quantum 
number, recent developments in optimized basis sets, particularly the use of atomic 
natural orbital contractions, have shown that results of full chemical accuracy (bet- 
ter than 1 kcal/mole even in the dissociation energy of N 2 ) can be achieved with 
manageable basis sets and reasonable MRCI expansions. 

In this review we have shown how FCI benchmark calculations can be used not 
only to provide general criteria by which approximate n-particle space treatments 
can be judged, but also to provide detailed calibration as to which treatments are 
appropriate in specific cases. In this way it is possible to attach confidence limits 
to the computed results obtained when such a treatment is applied in a large ANO 
basis. This ability to estimate realistic “error bars” for calculated quantities is very 
important when the results are to be used in interpreting experiment or in other 
calculations, and we can expect more use to be made of it in the future. We have 
described several applications of this technique to problems in spectroscopy (such 
as the D e values for NH and CN) and chemistry (F+H 2 barrier height). 

For the future, there is reason for considerable optimism about the method- 
ology and application of accurate quantum chemistry. A new generation of FCI 
programs should increase the size of FCI benchmark calculations by more than an 
order of magnitude, while the constant advances in more conventional methodolo- 
gies, even in areas thought to be exhausted, will continue. Further, new methods, 
particularly the application of multireference coupled-cluster approaches, can be ex- 
pected to find wider use, and some of the more exotic techniques that are currently 
being explored may begin to contribute to meaningful chemical applications. When 
combined with the wider availability and increasing power of supercomputers and 
minisupercomputers, the prospects for accurate quantum chemical calculations look 
bright indeed. 
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Table I. 1 Ai - 3 #i separation in CH 2 (kcal/mole) using a DZP basis set and corre 
lating the six valence electrons. 


Method 

Separation 

Error 

SCF° 

26.14 

14.17 

SCF“/SDCI 

14.63 

2.66 

SCF/SDCI+Q 

12.35 

0.38 

CPF 

12.42 

0.45 

tcscf 6 /mrci 

12.20 

0.23 

tcscf/mrci+q 

12.03 

0.06 

CASSCF C /MRCI 

11.97 

0.00 

CASSCF /MRCI+Q 

11.79 

-0.18 

FCI 

11.97 

— 


“ The SCF occupations are Iaj2aj3ail&2 

6 SCF treatment for Z B\ state, two-configuration MCSCF treatment for 1 A\ state 
(SCF configuration and 3 a\ — > 1 b\ excitation). 
c Active space comprises the C 2s 2p and H Is orbitals . 
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Table II. Spectroscopic constants for the 1 state of N 2 . 


Method 

SDCI 

SDCI+Q 

SDTCI 

SDTQCI 

SDQCI 

CPF 

MCPF 

MRCI 

MRCI+Q 

FCI 


r e ( 00 ) 

DZP basis 

6 electrons correlated 
w e (cm -1 ) 

D e (eV) 

2.102 

2436 

8.298 

2.115 

2373 

8.613 

2.107 

2411 

8.462 

2.121 

2343 

8.732 

2.116 

2361 

8.586 

2.112 

2382 

8.526 

2.114 

2370 

8.556 

2.123 

2334 

8.743 

2.123 

2333 

8.766 

2.123 

2333 

8.750 
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Table III. Polarizability of F ( al ) 


DZP + diffuse spd basis 
8 electrons correlated, a = d? E/dF 2 


Method 

a 

SDCI 

13.965 

SDCI+Q 

15.540 

CPF 

16.050 

MRCI® 

16.134 

MRCI+Q 

16.346 

MRCI'* 

16.034 

MRCI'+Q 

16.303 

FCI 

16.295 


“ The MRCI treatment is a CAS reference Cl based on a CASSCF calculation that 
included the 2 p electrons and the 2 p and 3 p orbitals in the active space. 
b The MRCI treatment is a CAS reference Cl based on a CASSCF calculation that 
included the 2s and 2 p electrons and the 2s, 2p, 3s and 3p orbitals in the active 
space. 
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Table IV. N( 4 5)/N 2 ( 1 S+ ) extended basis total energies (Eh), correlation ener- 
gies (E corr )(E h) and “dissociation energies” (eV). 


N atom 

Basis set 

Escf + 54 . 

E cor r 


(13s 8 p 6d) 

-0.400790 

-0.111493 


[6s 5 p 4 d] 

-0.400779 

-0.111321 


[5s 4 p 3d] 

-0.400769 

-0.110925 


[4s 3p 2d] 

-0.400725 

-0.109066 


(13s 8 p 6d 4/) 

-0.400790 

-0.121385 


[5s 4 p 3d 2/] 

-0.400769 

-0.120499 


[4s 3 p 2d 1/] 

-0.400725 

-0.117584 


[5s 4 p 3d 2/] (2 g) a 

-0.400769 

-0.122472 


[5s 4 p 3d 2/ 1^] 

-0.400769 

-0.122138 


N 2 molecule 6 

Basis set 

Escf + 108. E corr 

D e (SCF) D e (SDCI) 

(13s 8p 6d) 

-0.986307 -0.338118 

5.03 

8.16 

[6s 5 p 4 d] 

-0.985913 -0.337304 

5.02 

8.14 

[5s 4p 3d] 

-0.984833 -0.335395 

4.99 

8.08 

[4s 3p 2d] 

-0.983483 -0.329330 

4.95 

7.98 

(13s 8 p 6d 4/) 

-0.989318 -0.365735 

5.11 

8.45 

[5s 4p 3d 2/] 

-0.988031 -0.362548 

5.07 

8.38 

[4s 3 p 2d 1/j 

-0.986230 -0.353283 

5.03 

8.24 

[5s 4p 3d 2 j ](2g) a 

-0.988458 -0.370808 

5.09 

8.51 

[5s 4 p 3d 2/ 1#] 

-0.988322 -0.369270 

5.08 

8.48 


a 2 uncontracted g sets. 


6 r(N-N) = 2.1 do, 10 electrons correlated. 
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Table V. Theoretical study of the Mi - 3 Bi separation in CH 2 and SiH 2 . 


ANO basis set 

Separation (kcal/mole) a 

[3s 2 p Id /2s lp] 

ch 2 

11.33 

[4s 3 p 2d l//3s 2 p Id] 


9.66 

[5s 4p 3d 2/ 1<7/4s 3p 2d] 


9.24 

Expt+theory (T e ) 


9.28 (±0.1) 6 

Core- Valence 


+0.35 c 

Born-Oppenheimer 


- 0 . 1 41 

Relativistic 


-0.06“ 

[6s 5 p ld/3s lp] 

SiH 2 

-19.46 

[6s 5p 2d l//3s 2p Id] 


-20.15 

[6s 5p 3d 2/ lg/4s 3p 2d] 


-20.29 

Relativistic 


-0.30 

Zero- point 


— 0.33 e 

Expt(To) 


-21.0±0.7 / 


“ Ref. 88. 
6 Ref. 90. 
c Ref. 66. 
d Ref. 91. 
e Ref. 95. 
f Ref. 96. 
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Table VI. Summary of the Do values for the A 3 E state of NH, in eV. 



CH 

NH 

OH 

MRCP 

Recommended 0 

3.433 

3.344 

3.37±0.03 

4.360 

Experiment 6 

3.465 

>3.29 

<3.47 

4.393 


a Ref. 60. 

6 Refs. 108, 124 and 125. 


Table VII. Spectroscopic constants for the state of N 2 - 



[5s 

4 p 3 d 2/ 1^] basis 


Method 

r e (A) 

u> e (cm -1 ) 

D e (eV) 

MRCI(6) 

1.096 

2382 

10.015 

MRCI(6)+Q 

1.096 

2382 

10.042 

MRCI(10) 

1.101 

2343 

9.723 

MRCI(10)+Q 

1.102 

2336 

9.745 

Expt° 

1.098 

2359 

9.905 


“ Ref. 108. 
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Table VIII. Band positions and relative intensities of the HIR system of N2. 


Band 

Band positions (nm) 

Relative intensities 


Theory 

Experiment 0 

Theory 

Experiment 

0-0 

806 6 

806 

3.5 

4.0 

0-1 

854 

855 

1.0 C 

1.0 

0-2 

906 

907 

0.6 


0-3 

960 

963 

0.3 


0-4 

1016 

1023 

0.1 


1-0 

753 

752 

6.0 C 

6.0 

1-1 

795 

795 

0.2 


1-2 

840 

840 

0.7 

<0.5 

1-3 

885 

887 

1.6 


1-4 

932 

938 

1.4 


2-0 

707 

706 

I — 1 

O 

O 

O 

10.0 

2-1 

745 

744 

8.8 

8.0 

2-2 

783 

783 

6.6 

5.0 

2-3 

823 

824 

0.3 

<0.5 

2-4 

864 

868 

0.7 


3-0 

668 

667 

1.5 C 

1.5 

3-1 

701 

700 

10.4 

9.0 

2-3 

735 

735 

0.6 

<0.5 

3-3 

771 

771 

4.8 

3.0 

3-4 

806 

809 

2.5 



“Ref. 152. 

fc The 0-0 bands were shifted into coincidence; this required a shift of 567 cm 1 in 
the theoretical T e . 

c These bands were normalized to experiment by adjusting the vibrational popula- 
tions. This requires that the relative populations of v'=0-3 be 0.073, 0.295, 1.00, 
and 0.725, respectively. 
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Table IX. FCI calibration of the classical barrier height of F+H 2 — >HF+H a . 


A. At the FCI saddle point 


FCI 

MRCI(300) 

MRCI(300)+Q 

MRCI(322)(0.05) 

MRCI(322)(0.05)+Q 

MRCI(322)(0.025) 6 

MRCI(322)(0.025)+Q 6 

ACPF(322)(0.025)(12-Ref) 6 

ACPF(322)(0.025)(13-Ref) c 

MRCI(322)(0.01) 

MRCI(322)(0.01)+Q 

MRCI(522)(0.025) 

MRCI(522)(0.025)+Q 


barrier 

exothermicity 

4.50 

28.84 

5.18 

28.57 

4.43 

29.12 

5.00 

29.12 

4.32 

29.21 

4.73 

29.17 

4.51 

28.80 

4.56 

28.89 

4.57 

28.89 

4.71 

29.19 

4.54 

28.84 

4.55 

29.41 

4.32 

29.31 


B. At the optimized saddle-point geometry d 


FCI 

MRCI(300) 

CPF 

MRCI(300)+Q 

MRCI(322)(0.025) 6 

MRCI(322)(0.025)+Q 6 

ACPF(322)(0.025) C 

SDTQCI e 


(F-H) 

r(H-H) 

AE b 

AE rx 

2.761 

1.467 

4.50 

28.84 

2.740 

1.476 

5.16 

28.57 

2.801 

1.467 

4.40 

26.47 

2.795 

1.467 

4.42 

29.12 

2.761 

1.474 

4.70 

29.17 

2.755 

1.475 

4.49 

28.80 

2.760 

1.475 

4.56 

28.89 

2.763 

1.465 

4.45 



“ Energies in kcal/mole and bond lengths in a 0 . All calculations are done using the 
[4s 3 p Id /2s 1 p] basis set and correlating seven electrons. The barrier is referenced 
to F. . .H 2 (50<zo), and the exothermicity is computed using HF. . .H(50ao)- 
6 The 12 reference configurations chosen from the CASSCF wave function. 
c Includes the additional configuration found to be important in the ACPF calcu- 
lation — see the text. 

^ Geometry optimized using a biquadratic fit to a grid of nine points. 
e Scuseria and Schaefer, Ref. 173. 
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Table X. Theoretical studies of the classical saddle-point geometry (a 0 ) and ener- 
getics (kcal/mole) for the F+H 2 reaction. 


Basis® 

12-Ref d 

A 

A 

A 

A 

A 

A 

A 

A 


A 

A 

a+h (fy 

A+H (/) 
A-F (g) 
A— F(sf) 

13- Ref' 1 
A 
A 
A 
A 
A 
A 


Level of treatment 


MRCI(322)(2p) e 

MRCI(322)(2p)+Q 

MRCI(322)(2p) 

MRCI(322)(2p)+Q 

ACPF(322)(2p) 

MRCI(322) 

MRCI(322)+Q 

ACPF(322) 

CCI(322) 

CCI(322)+Q 

CCI(322) 

CCI(322) +Q 

CCI(322) 

CCI(322)+Q 


MRCI(322)(2p) 

MRCI(322)(2p)+Q 

ACPF(322)(2p) 

MRCI(322) 

MRCI(322)+Q 

ACPF(322) 


r 6 

r HF 

-6 

r HH 

baxrier c 

exothermicity 

2.899 

1.455 

2.99 

33.96 

2.910 

1.456 

2.42 

33.42 

(2.95) 

(1.45) 

2.96 


(2.95) 

(1.45) 

2.40 


(2.95) 

(1.45) 

2.45 


2.914 

1.451 

2.63 

31.61 

2.950 

1.450 

1.66 

30.47 

(2.95) 

(1.45) 

1.17 



.../ .../ 2.79 31.8 

.../ ...f 2.02 30.7 

.../ .. / 2.73 


J .J 1.95 


2.879 

1.447 

2.89 


2.909 

1.445 

2.14 


(2.95) 

( L45 ) 

2.96 

33.96 

2.95) 

(1.45) 

2.51 

33.42 

2.914 

1.453 

2.61 

33.54 

(2.95) 

(1.45) 

2.59 

31.61 

2.95) 

(1.45) 

1.81 

30.47 

2.967 

1.447 

1.85 

30.58 


Expt.' 


31.73 


“ “A” denotes the [5s 5 p 3 d 2/ lg/4s 3 p 2d] basis described in Ref. 168. 

6 The saddle-point geometries in parentheses have not been optimized. 
c The barrier is referenced to F. . .H2(50ao)i nnd the exothermicity is computed 
using HF. . .H(50ao). 

d The 12 reference configurations chosen from the CASSCF wave function. 
e (2p) indicates a seven-electron treatment (i.e. 2s correlation is excluded). 

/ The MRCI(300)+Q saddle point geometry is used, r(F-H)=2.921 a 0 and r(H- 
H)=1.450 a 0 . 

9 Denotes that a function of this angular momentum type has been added. 

^ Includes the additional configuration found to be important in the ACPF calcu- 
lation — see the text. 

1 From data in Ref. 108. 
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Table XI. Study of the 2,< effect on the F+H 2 barrier height(kcal/mole)“ in the 
DZ +pf basis set. 


Level of treatment 

7 electron 
9 electron 
A 


MRCI 

MRCI+Q 

ACPF 

SDTQCI 

FCI 

8.73 

8.74 

8.75 

8.82 

8.77 

6.63 

6.50 

6.53 

6.82 

6.64 

2.10 

2.24 

2.22 

2.00 

2.13 


° The saddle point geometries are taken from the MRCI(2p)+Q and MRCI+Q 
calculations in this basis set. 
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Table XII. Zero-point and tunnelling effects on the barrier height of the F + H 2 
and F + D 2 reactions. 


F + H 2 surface 


a o 

r HHi a 0 

Barrier, kcal/mole 
Sym. stretch, cm -1 
Bend, cm -1 
Asym. stretch®, cm 1 
Zero-point correction 6 , kcal/mole 
E barrier + zero point, kcal/mole 
Tunnelling correction, kcal/mole 
Threshold, kcal/mole 


Classical Barrier 


CCI 

CCI + Q 

2.879 

2.909 

1.447 

1.445 

2.888 

2.143 

3706 

3768 

68.5 

45.9 

692* 

605* 

-0.602 

—0.643 

2.286 

1.500 

-0.54 

-0.47 

1.75 

1.03 


Adiabatic Barrier 


CCI 

CCI + Q 

3.070 

3.155 

1.425 

1.421 

2.639 

1.860 

4074 

4178 

68.5 

45.9 

530* 

371* 

-0.076 

-0.057 

2.563 

1.803 

-0.42 

-0.29 

2.14 

1.51 


F + D 2 surface 


THF, Oo 

Barrier, kcal/mole 
Sym. stretch, cm -1 
Bend, cm -1 
Asym. stretch®, cm -1 
Zero-point correction 6 , kcal/mole 
E barrier + zero point, kcal/mole 
Tunnelling correction, kcal/mole 
Threshold, kcal/mole 


Classical Barrier 


CCI 

CCI -1- Q 

2.879 

2.909 

1.447 

1.445 

2.888 

2.143 

2623 

2667 

37.7 

19.1 

512z 

448* 

— 0.488 

-0.532 

2.400 

1.611 

-0.40 

-0.35 

2.00 

1.26 


Adiabatic Barrier 


CCI 

CCI + Q 

3.010 

3.075 

1.430 

1.427 

2.761 

1.997 

2811 

2876 

37.7 

19.1 

428*' 

334* 

-0.220 

-0.233 

2.541 

1.764 

-0.34 

-0.26 

2.20 

1.50 


“From the normal mode analysis at the classical barrier, and computed from the 
curvature along the Eckait potential at the adiabatic barrier. 

6 For H 2 (D 2 ) we used u; e = 4401(3116) cm -1 , respectively, from Ref. 108. 
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Figure Captions 

Figure 1: Comparison of the theoretical radiative lifetimes (63) for different ro- 
tational levels of the u=0 level of the A 2 Yj+ state of OH (dashed line) with the 
high frequency deflection measurements (F x component) of Brzozowski et at. (116) 
shown with error bars. 

Figure 2. The FCI potential energy curves for the X 1 E + , A 1 II, and C E + states 
of A1H (68). 

Figure 3. Potential energy curves and vibrational levels for the A 3 E+ , A' 5 E+ and 
B 3 U g states of N 2 (150). 
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